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We study post weak decoupling coherent active-sterile and active-active matter-enhanced neutrino 
flavor transformation in the early universe. We show that flavor conversion efficiency at Mikheyev- 
Smirnov-Wolfenstein resonances is likely to be high (adiabatic evolution) for relevant neutrino pa- 
rameters and energies. However, we point out that these resonances cannot sweep smoothly and 
continuously with the expansion of the universe. We show how neutrino flavor conversion in this 
way can leave both the active and sterile neutrinos with non-thermal energy spectra, and how, in 
turn, these distorted energy spectra can affect the neutron-to-proton ratio, primordial nucleosyn- 
thesis, and cosmological mass/closure constraints on sterile neutrinos. We demonstrate that the 
existence of a light sterile neutrino which mixes with active neutrinos can change fundamentally the 
relationship between the cosmological lepton numbers and the primordial nucleosynthesis 4 He yield. 

PACS numbers: 14.60.Pq; 14.60.St; 26.35,+c; 95.30.-k 



I. INTRODUCTION 

If light sterile neutrinos exist we will be forced to 
re-think the role of the weak interaction in the early 
universe, primordial nucleosynthesis, and cosmology. 
Though light sterile neutrinos which mix with active neu- 
trinos long have been a subject of theoretical speculation, 
the LSND and KARMEN experiments Q gave rise to 
particular interest in the mass-squared difference range 
0.2 eV 2 < Sm 2 s < 100 eV 2 . Here we study for this range 
of <5m 2 s the cosmological lepton number-driven conver- 
sion of active neutrinos, v a (and/or u a ) with a = e, y,, r, 
to a singlet, "sterile" neutrino species v s (or D s ) in the 
epoch of the early universe after decoupling of the weak 
interactions, when neutrino spectral distortions are likely 
to persist. 

A positive signal in the on-going mini-BooNE experi- 
ment 0, «.e., confirming the interpretation of the LSND 
result in terms of vacuum neutrino mixing, sets up an im- 
mediate crisis in neutrino physics. Such a result, when 
combined with the already well established evidence for 
neutrino mixing at mass-squared differences associated 
with the atmospheric (5m 2 ~ 3 x 10~ 3 eV 2 ) and solar 
neutrino (5m 2 ~ 7x 10~ 5 eV 2 ) anomalies, would suggest 
the existence of three independent neutrino mass-squared 
differences which would, in turn, require four neutrino 
species. Given the Z°-width limit on the number of fla- 
vors of neutrinos with standard weak interactions (3), 
a fourth neutrino would have to be "sterile," with sub- 
weak interaction strength, e.g., perhaps an SU(2) singlet. 
The only alternative to this line of reasoning and to this 



conclusion is the possibility of CPT violation |3|. How- 
ever, there is no consistency of the neutrino oscillation 
data with a CPT-violating three-neutrino model at a 3-cr 
level Q. 

Hand in hand with this particle physics dilemma, evi- 
dence for a singlet neutrino that mixes with active neutri- 
nos in this mass-squared range also confronts cosmology 
with a curious and vexing problem. In the standard cos- 
mological model with zero or near-zero net lepton num- 
bers one would expect that matter-supressed neutrino 
oscillations in the channel v a ^ v s or in D a ^ D s (where 
a = e, /i, r) proceeding in the regime above weak interac- 
tion decoupling (T > 3 MeV) , would efficiently populate 
seas of singlet neutrinos |£g. The significant additional 
energy density in these sterile neutrino seas would engen- 
der a faster expansion rate for the universe and a conse- 
quently higher temperature for Weak Freeze-Out (where 
the initial isospin of the universe, the neutron-to-proton 
ratio is set). A higher Weak Freeze-Out temperature 
would result in more neutrons and, hence, a higher yield 
of 4 He. 

A higher predicted abundance of 4 He arguably may be 
in conflict (or is close to being in conflict) with the obser- 
vationally inferred upper limit on the primordial helium 
abundance. Depending on the helium abundance inferred 
from compact blue galaxies, an increase in the predicted 
Big Bang Nucleosynthesis (BBN) 4 He yield may or may 
not be disfavored 1^, . However, the primordial helium 
abundance is notoriously difficult to extract from the ob- 
servational data and recent studies point to a fair range 
for the observationally inferred primordial mass helium 
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fraction: 23% to 26% 7]. The upper limit of this range is 
provocatively close to the standard BBN 4 He mass frac- 
tion yield prediction, 24.85 ± 0.05%, as computed with 
the deuterium-determined or CMB (Cosmic Microwave 
Background) anisotropy-determined baryon density. 

Additionally, it has been suggested that a fully pop- 
ulated sea of sterile neutrinos and antineutrinos with rest 
masses ~ y <5m 2 s could be in conflict with neutrino mass 
bounds derived from CMB anisotropy limits and large 
scale structure considerations . There is a recent analy- 
sis of constraints from measurements of galaxy bias stem- 
ming from galaxy-galaxy lensing and the inferred linear 
matter power spectrum derived from the Lyman alpha 
forest in the Sloan Digital Sky Survey (SDSS) p]. This 
analysis specifically considers a so-called "3+1" neutrino 
mass hierarchy, i.e., the scheme which is appropriate for 
constraining sterile neutrinos. The neutrino mass con- 
straint so derived is somewhat less stringent than con- 
straints in schemes with three neutrinos with degenerate 
masses. However, the central conclusions of Ref. |8j sur- 
vive. 

Should we someday be confronted with a positive indi- 
cation of neutrino flavor mixing with mass-squared scale 
consistent with the range for <5to 2 s , we will have a prob- 
lem that would call for modification either of our notions 
of basic neutrino physics or of the standard cosmological 
model. There have been a number of ways proposed to 
get out of these cosmological difficulties. For instance, if 
neutrinos do not acquire mass until after the BBN epoch, 
as may occur via a late-time phase transition pTlll2j| . the 
singlet states will not be populated via oscillations dur- 
ing the BBN era. In addition to schemes involving the 
epoch of neutrino mass generation |l3j . annihilation 
or decay [l2T | of the singlet neutrinos (when m s ~ T) may 
alleviate or avoid the CMB and large scale structure con- 
straints. However, chief among the mechanisms proposed 
to escape the cosmological difficulties associated with sin- 
glet neutrinos, is the the invocation of a significant net 
lepton number [l5|. 

The idea is that the net lepton number gives active 
neutrinos larger effective masses in medium in the early 
universe, thereby driving them further off-resonance in 
the epoch prior to weak decoupling (i.e., T > 3 MeV) and 
reducing their effective matter mixing angles with the 
singlet neutrino. In turn, smaller effective matter mixing 
angles would imply a suppressed production of singlet 
neutrinos and, hence, a reduced population of the singlet 
neutrino sea. This lepton number-induced suppression 
of active-sterile mixing at high temperature is why we 
assume here that there is no initial population of the 
sterile neutrino sea. 

The lepton number residing in the sea of v a and v a 
neutrinos (a = e, /i, r) is defined in analogy to the baryon 
number rj = [rib — ng)/n 7 , 



n Va - np Q 



(1) 



density at temperature T 7 , and where n Va and n„ a are 
the number densities of v a and v a neutrinos, respectively, 
at this epoch. After the epoch of e ± annihilation the 
baryon number is r\ « 6x 10 -10 , whereas at earlier epochs 
it is roughly two and half times larger. We will consider 
here net lepton numbers which are vastly larger than 77, 
so its precise value is of no consequence for our results. 

We can insure that the effective matter mixing angles 
for the oscillation channel v a # v s (or v a ^ v s ) are 
sufficiently small to suppress singlet neutrino produc- 
tion if the Mikheyev-Smirnov-Wolfenstein (MSW) pij 
resonance temperature is less than the weak decoupling 
temperature, T res < 7d ec . This implies that the lepton 
number associated with any of the active neutrino flavors 
should satisfy, 
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where n 7 = (2C(3)/tt 2 ) T 3 is the proper photon number 



where 9 is the vacuum mixing angle characteristic of 
v a v a oscillations, iVdegen is the number of neu- 
trino species possessing this lepton number, and where 
e = E v /T. For neutrinos with typical energies in the 
early universe (i.e., e ~ 1), suppression of singlet neu- 
trino production would require lepton numbers ranging 
from L > 10-* for 5m 2 as = 0.2 eV 2 to L > 5 x 10~ 3 for 
Strips = 10 eV 2 . Current limits on lepton numbers are 
\L Ua I < 0.1 [17J (and possibly even weaker by a factor of 
two or so if allowance is made for another source of extra 
energy density in the early universe [l8L[T9l |'). Therefore, 
this avenue for escape from the sterile neutrino conun- 
drum appears to be allowed, albeit at the cost of a huge 
disparity between the lepton and baryon numbers. 

However, this argument overlooks an important point. 
Though the large lepton number supresses the effective 
matter mixing angle for v a ^ v s during the epoch 
of the early universe where active neutrinos are ther- 
mally coupled (T > Tdcc), it can cause coherent matter- 
enhancement of this channel at lower temperatures where 
the active neutrinos rarely scatter and are effectively de- 
coupled. Resonant MSW transformation of active neu- 
trinos to singlets in the channel v a ^ v a is, however, 
self limiting. This is because as the universe expands 
and the resonance sweeps from low toward higher neu- 
trino energy, the conversion of v^s decreases the lep- 
ton number which, in turn, causes the resonance sweep 
rate to increase, eventually causing neutrinos to evolve 
non-adiabatically through resonance and so causing fla- 
vor transformation to cease. 

At issue then is how many active neutrinos can be con- 
verted to sterile neutrinos prior to or during the epoch 
where the neutron-to-proton ratio is set ( "Weak Freeze- 
Out"). If there is a significant conversion, the resultant 
non-thermal active neutrino energy spectra can cause an 
increase or decrease (if v a v s is enhanced) in the 4 He 
yield and call into question the viability of invoking a 
large net lepton number to reconcile neutrino physics 
and BBN. Other but related aspects of transformation- 
induced nonthermal neutrino spectra effects on primor- 



3 



dial nucleosynthesis have been studied in Ref.s pol. l2l|. 
Likewise, previous studies have considered other aspects 
of the relationship between sterile neutrinos and BBN 
|22|. as well as constraints on sterile neutrinos without a 
primordial lepton number |23j . In any case, non-thermal 
energy distribution functions for v e and/or v e change the 
relationship between the BBN 4 He yield and the neutrino 
chemical potentials. 

In section II we discuss the physics of active-sterile 
neutrino flavor transformation in the early universe and 
point out a key issue in how the MSW resonance sweeps 
through the neutrino energy distribution functions as the 
universe expands. The generally high adiabaticity of neu- 
trino flavor evolution is also pointed out in this section. 
Simultaneous active-active and active-sterile neutrino fla- 
vor conversion, and "synchronization" are also discussed 
in this section. Possible multi-neutrino mass level cross- 
ing scenarios in the early universe are discussed in this 
section. Sterile neutrino contributions to closure, con- 
straints on this from large scale structure and Cosmic 
Microwave Background radiation considerations, as well 
as other sterile neutrino sea population constraints are 
examined in section III. In section IV we describe how 
distorted v e and/or v e distribution functions impact the 
rates of the lepton capture reactions that determine the 
neutron-to-proton ratio and the 4 He yield in BBN. This 
is then applied in various initial lepton number and neu- 
trino conversion scenarios. Finally, in section V we give 
conclusions and speculations regarding the neutrino mass 
and cosmological lepton number insights that would fol- 
low in the wake of an experimental signature for a large 
neutrino mass-squared difference of order the range given 
for 5m^ s . Appendix A provides an exposition of the lep- 
ton capture rates on free nucleons when, as appropriate, 
v e or D e energy distribution functions are zero up to some 
energy, and thermal/Fermi-Dirac at higher energies. 

II. COHERENT NEUTRINO FLAVOR 
TRANSFORMATION IN THE EARLY UNIVERSE 

Coherent conversion of active neutrino species into sin- 
glets in the early universe can occur through the usual 



Here n e = n e - — n e + is the net number density of elec- 
trons, n n = rib — Tip is the number density of neutrons, 



MSW process, albeit in an exotic setting. This pro- 
cess can be described simply: (1) an active neutrino 
(mostly the light mass state in vacuum) forward scat- 
ters on particles in the plasma and, if there is a net 
lepton and/or baryon number, will acquire a positive ef- 
fective mass; (2) an MSW resonance (mass level cross- 
ing, where in-medium mixing is large) can occur when 
this effective mass is close to the mass associated with 
the singlet (mostly the heavy mass state). This process 
is discussed in subsection A, but despite the simplic- 
ity of the physics behind it, the neutrino energy depen- 
dence/history of MSW resonances in the early universe 
can be quite complex, as shown in subsection B. 

The efficiency of flavor conversion at a mass level cross- 
ing depends on the ratio of the resonance width (in time 
or space) to the neutrino oscillation length. Efficient, 
adiabatic conversion takes place only when this ratio is 
large, and in subsection C we examine this physics in 
detail for our particular problem. 

Active-active matter-enhanced neutrino flavor conver- 
sion, discussed in subsection D, can occur simultaneously 
with active-sterile transformation in the early universe. 
This can greatly complicate computing the history of the 
neutrino distribution functions. Subsection E below deals 
with the limit where active- active conversion is efficient, 
while subsection F examines flavor evolution in the limit 
where this conversion channel is inefficient. 



A. Neutrino Effective Masses and Level Crossings 

The forward charged and neutral current exchange 
Hamiltonians for the neutrinos in the early universe are 
as follows (see, e.g., |23): 

J3>,) = (3) 



(4) 



(5) 



(6) 



and rib and n p are the net number densities of baryons 



H (v e ) = V2G F \n e - -n n J + s/2G F (2 (n„ e - n p J + (n V)l - n p J + {n Vr - n Pr )) 
H [vp) = V2G F (~\ n nj + V2G F ((n„ e - n p J + 2 (n Vft - n p J + (u Vt - m>J)) 
H{v T ) = V2G F (~-n n j +V2G F ((n Ve - n p J + (rc^ - n P J + 2 (n„ r -n p J). 
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and protons, respectively. Charge neutrality implies that 
the number density of protons is n p = n e = ni,Y e . The 
net number of electrons per baryon is Y e . The baryon 
number density is n.& rs r/n^, where the baryon-to-photon 
ratio rj is as defined above. 

Weak Decoupling occurs when neutrino scattering be- 
comes so slow that it can no longer facilitate efficient en- 
ergy exchange between the neutrino gas and the plasma. 
For the low lepton numbers considered here, Weak De- 
coupling occurs around temperature T ~ 3 MeV, though 
this decoupling process takes place over a range in tem- 
perature of a few MeV. 

Weak Freeze- Out occurs when the rates of the reac- 
tions that govern the ratio of neutrons-to-protons (n/p = 
1/Y e — 1) fall below the expansion rate of the universe. 
This is usually taken to be T w 0.7 MeV for standard cos- 
mological parameters. However, below this temperature 
Y e continues to be modified by lepton capture and/or free 
neutron decay as discussed below. 

Note that for temperatures well above Weak Decou- 
pling, we have Y e w 0.5, i.e., nearly equal numbers 
of neutrons and protons. We could then approximate 
n e -\n n = nh(§V e -±) w n 7 ?//4, and -|n„ = n b (Y e /2- 
i) Rj — n 7 7//4. We will use this approximation in what 
follows even for the epoch below weak decoupling where 
it is not numerically accurate. This will result in no loss 
of accuracy in the full calculation because we consider 
large net lepton numbers, L Va ^> r\. 

We can denote the weak potentials from neutrino- 
electron charged current forward exchange scattering and 
neutrino-neutrino neutral current forward exchange scat- 
tering as A and B, respectively, with their sum being 



A + B 



2V2C(3)G F T 3 



(7) 



We will refer to this quantity as the "potential lepton 
number." In general this may be different for different 
channels v a v s , even for a given set of lepton numbers 
associated with each flavor. 

Finally, since the early universe is at relatively high en- 
tropy per baryon, the overall weak potential has a contri- 
bution from neutrino neutral current forward scattering 
on a thermal lepton background. This thermal potential 



C « -r a G 2 F eT\ (9) 
where the neutrino energy divided by the temperature 
is e = E„/T. For the conversion channel v e v s , we 



employ 
use r,, ^ 



79.34, while for the channel v 



22.22. If the neutrinos have strictly thermal 



energy distribution functions, then 



F 2 (0) F 2 (0) 



(10) 



where the neutrino and antineutrino degeneracy parame- 
ters are rj Va and r\v a , respectively, and the Fermi integrals 
of order 2 arc defined below. The approximations r e rs r e 
and r M , r » r? r suffice for lepton numbers below the con- 
ventional limits. 



The total weak forward scattering potential is 



V = A + B + C. 



(11) 



For the transformation channel v a v s , the neutrino 
mass level crossing (MSW resonance) condition for a neu- 
trino with scaled energy e is 



where Gf is the Fermi constant, the Riemann Zeta func- 
tion of argument 3 is ( (3) « 1.20206, and we take the 
plus sign for transformation of v e , and the minus sign for 
conversion of and/or v T . (Here the plus sign is taken 
when we intend A + B = H(u e ) and the minus sign is 
taken when A + B = H{v^ iT ).) A measure of the lepton 
number which enters into the potential for the v a ^ v s 
(a = e, (i. t) conversion channel is 



C = 2L^ 



(8) 



Sm 2 cos 29 

2lr 



v, 



(12) 



where 5m 2 is the difference of the squares of the appro- 
priate neutrino mass eigenvalues and 9 is the relevant 
effective two-by-two vacuum mixing angle. Neglecting 
the light mass eigenvalue, the effective mass-squared ac- 
quired by an electron neutrino from forward scattering 
on weak charge-carrying targets in the early universe is 



"eff 



i 2eV« (8.03 x 10 



MeV 2 ) e(C± V /A)(^j - (2.16 x 



HT 20 MeV 2 ) e 2 ( —\ 



(13) 



It is clear that we can negelect the second term (the ther- 0.7 MeV. We also neglect the baryon/electron term, 
mal term C) in Eq. (|llfl in the regime between Weak ±^/4. 
Decoupling and Weak Freeze-Out, where 3 MeV > T > 
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At a given temperature, the scaled neutrino energy 
which is resonant is then 



5m 2 cos 29 
2VT 



(14) 



The dependence of resonant neutrino energy on temper- 
ature and lepton number is 



ir 2 6m 2 cos 29 



2 5 /2^(3) Gf(£±?7/4)T 4 
0.124| j — j — i 



(15) 



leW 



C 



It is clear from Eq. (|15l) that as the universe expands 
and the temperature drops, the resonance energy e rcs will 
sweep from lower to higher values. In fact, as the reso- 
nance sweeps through the active neutrino distribution, 
converting v a — > v s , C will decrease, further accelerating 
the resonance sweep rate. 

Assuming homogeneity and isotropy, the number den- 
sity of active neutrinos v a with thermal distribution func- 
tion f Va (e) in the scaled energy range e to e + de is 



(16) 



where n Va is the total number density (that is, integrated 
over all neutrino energies). In terms of the temperature T 
and degeneracy parameter rj Ua = [i Vct /T, where \i Va is the 
appropriate chemical potential, the thermal distribution 
function is 



fu a (e) = 



1 



F 2 {ilu a ) e'-^c +1 



(17) 



We define relativistic Fermi integrals of order k in the 
usual fashion: 



F k fa) = 



v h dx 



e x- v + i 



(18) 



The total number density of thermally distributed ac- 
tive neutrinos v a with temperature T v and degeneracy 
parameter rj Va is 



rp3 
_ V 

2tt 



(19) 



Note that if the neutrino degeneracy parameter is r\ Va — 
0, then F 2 (0) = 3£ (3) /2 and the number density of ther- 
mally distributed v^s is 



(20) 



where we allow for the neutrino temperature T v to differ 
from the photon/plasma temperature T 1 . 

The relationship between the lepton number in a flavor 
neutrinos and the v a degeneracy parameter is 



(wk)(^) 3 ["- + <" 2 ]- (21) 



This relation assumes that neutrinos have Fermi-Dirac 
energy spectra and that v a and antineutrinos v a are (or 
were at one point) in thermal and chemical equilibrium so 
that r\ij a = —r\v a ■ In the limit where the lepton number 
is small, so that r] Va <^ 1, and the neutrino and photon 
temperatures are nearly the same, we can approximate 
Eq. I|21|) as 7y„ Q sa lAQL Wa . Neutrino degeracy parameter 
is a comoving invariant; whereas, lepton number is not 
in general since the photons can be heated relative to the 
neutrinos by, e.g., annihilation. 



B. Lepton Number Depletion and the 
Time/Temperature Dependence of Resonance 
Energies 

As the universe expands and v a neutrinos are con- 
verted to sterile species v Sl the lepton number L Va drops. 
As C approaches zero, the resonance sweep rate becomes 
so large that neutrinos will be propagating through MSW 
resonances non-adiabatically 25]. Efficient neutrino fla- 
vor conversion ceases at this point. If the conversion 
process results in a change in the number density of v a 
neutrinos, An^ , such that the lepton number associated 
with this species changes by AL„ Q = ~IS.n Va jn 1 , then 
the potential lepton number would change from its ini- 



tial value, C 



initial 



to 



^final £\ 



titial 



2AL„ 



(22) 



The adiabaticity condition ensures that flavor conversion 
ceases when £ final approaches zero, or in other words, 
when AL Un = — £ lmtlal /2. It is important to note that 
transformation of any flavor active neutrino to sterile fla- 
vor can drive down the overall potential lepton number, 
no matter which flavor or flavors of active neutrinos har- 
bor the net lepton number. 

If additionally we were to assume that the resonance 
smoothly and continuously swept through the v a energy 
distribution from zero to scaled energy e during this con- 
version process, we would have An„ a = L dn Ua and so 
the concomitant change in lepton number would be 



Try 



1 



Fi (0) Jo e*-^ a + 1 



-dx. 



(23) 



In this idealized scenario, the potential lepton number as 
a function of e is 



initial 



3 l 



4 F 2 (0) J e x ~^ a + 1 



-dx. 



(24) 



In this last relation we have set the photon/plasma and 
neutrino temperatures to be the same. This is a good 
approximation in the epoch where it turns out we will be 
most interested in resonance sweep, between Weak De- 
coupling and T « 0.5 MeV. During this time there has 
been little annihilation of e pairs and, consequently, lit- 
tle heating of the photons/plasma relative to the decou- 
pled neutrinos. 
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FIG. 1: The nonthermal scaled energy (E„/T) distributions 
/ (E v /T) for v s (dashed) and v £ (solid) resulting from smooth, 
adiabatic resonance sweep from E v jT = to E v /T = e. 

Employing the approximation of a smooth and contin- 
uous sweep of scaled resonance energy from zero to e, we 
can re- write the resonance condition, Eq. (|15|l . as 

rf \~ n 2 6m 2 cos28 
££(e) * 2 5 / 2 C (3) G-fT a ' (25) 

For smooth, continuous and adiabatic (i.e., complete 
conversion) resonance sweep up to a scaled energy e, the 
resulting active and sterile neutrino distribution func- 
tions would be as shown in Fig. JQ). This energy dis- 
tribution function is zero for all values of scaled neutrino 
energy < E v /T < e, and has a normal Fermi-Dirac 
thermal distribution character for all neutrino energies 
E v /T > e. The corresponding sterile neutrino energy 
spectrum would be the "mirror image" of this: a ther- 
mal Fermi-Dirac spectrum for < E v /T < e, and zero 
occupation for E v jT > e. 

However, Eq. I|25fl reveals a problem: the resonance 
cannot sweep continuously and smoothly to the point 
where £(e) — > 0. This is because e£(e) is a peaked func- 
tion. The maximum of this function occurs for a value 
£max satisfying the integral equation 

4 ax » 2( (3) (e e °» + l)£(e max ). (26) 

It is clear, however, that as the universe expands, the 
right hand side of Eq. (|25|l will increase monotonically. 
Although the resonance sweep can begin smoothly and 
continuously, there will come a point where it is no longer 
possible to find a solution to Eq. Q25J1. This will occur 
when the resonance energy reaches e max . 

Fig. J2J shows graphically the problem of obtaining 
a solution to Eq. l|25[l for a particular case. The solid 
line in this figure is e£(e) as computed by assuming a 
smooth and continuous resonance sweep scenario. Here 
we have chosen initial lepton numbers L„ c = L u = 
L Vt = 0.01096. This corresponds to an initial potential 
lepton number £(e = 0) = 0.04384. The arrows in this 



figure give the sense of evolution along the solid curve as 
the universe expands and the net potential lepton num- 
ber decreases as a result of neutrino flavor conversion in 
the channel u e — > v s . The maximum value on this curve 
occurs at e max sa 0.598. 

However, the potential lepton number would be com- 
pletely depleted (i.e., £(e c o .) = 0) in the smooth and 
continuous resonance sweep scenario only when e reaches 
a "cut-off" value e c . ., which is 0.987 in this case. 
The horizontal dashed lines in Fig. @ correspond to val- 
ues of the right-hand side of Eq. (|25() for epochs of the 
universe corresponding to temperatures T — 2.0 MeV, 
T = 1.595 MeV, and T = 1.5 MeV. Solutions to this 
equation are possible when these curves cross the solid 
e-L(e) curve. Clearly, no solutions are possible in the con- 
tinuous resonance sweep scenario for T < 1.6 MeV in this 
case. 

What happens beyond this point, e.g., for T < 
1.6 MeV? If we relax the demand that the resonance 
sweep be continuous, then it is possible in principle to 
find a solution to Eq. I|25|l as the temperature drops be- 
yond the point where e = e max , though this would require 
that the product of scaled resonance energy and potential 
lepton number differ from the solid curve e£(e). Though 
a detailed numerical model of this process is beyond the 
scope of this work, we can get a rough idea of what might 
happen with the following argument. 

Suppose we take a time step resulting in a new temper- 
ature T" slightly lower than T max , the temperature where 
the last continuous sweep solution exists, i.e., where 
e = fmax (T max ~ 1.6 MeV for our example case) and 
£ = £(e max ). One possibility is that the resonance en- 
ergy could skip to some value e > e max , toward the higher 
energy portion of the neutrino distribution function. In 
this way the product e£(e ma x) could be large enough to 
match the right-hand side of Eq. I|25[) at the new tem- 
perature T". Of course, this would result in resonant 
neutrino flavor conversion and so £ would be lowered 
and eventually we again would be unable to maintain a 
smooth resonance sweep. At that point the resonance en- 
ergy could skip again discontinuously, etc.. It is possible 
that beyond e max the resonance sweeps stochastically in 
this way through relatively small intervals of the active 
neutrino distribution function leaving a "picket fence" 
distribution beyond e max . 

Though the details may differ from this simple scheme, 
we believe that resonant neutrino flavor conversion for 
e > fmax will occur because: (1) as we show below, neu- 
trino flavor evolution at this epoch is very adiabatic for 
relevant neutrino parameters; and (2) the resonance con- 
dition can be met for some value of neutrino energy so 
long as the net lepton number is non-zero. In any case, 
however, active-to-sterile neutrino conversion v a — > v s 
will have to cease when £ approaches zero. 

At this point we will be left with grossly non-thermal, 
non- Fermi-Dirac v a and v s distributions. Since this pro- 
cess occurs after Weak Decoupling, active neutrino in- 
clastic scattering processes on electrons, nucleons, and 
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FIG. 2: The solid line is eC(e), the product of e and potential 
lepton number, in the smooth and continuous resonance sweep 
case for initial lepton numbers L„ e = L v = L Vt — 0.01096, 
corresponding to initial potential lepton number C(e = 0) = 
0.04384. The arrows give the sense of evolution along this 
curve as the universe expands and the net potential lepton 
number decreases as a result of neutrino flavor conversion 
in the channel v e — > v s with 8m cos 26» = leV 2 . The hor- 
izontal dashed lines correspond to values of the right-hand 
side of Eq. (I25H for the indicated epochs (temperatures). So- 
lutions to Eq. i12-")l are possible at a given epoch when the 
corresponding dashed line crosses the e£(e) curve. Physical 
solutions are circled here for T = 2.0 Me V and 1.6 Me V. The 
maximum value on the e£(e) curve occurs at e ma x ~ 0.598, 
and in the smooth resonance sweep scenario this is reached 
at T w 1.6 MeV. Clearly, no solutions are possible in this 
scenario for T < 1.6 MeV. If the system were forced to follow 
the smooth resonance sweep e£(e) curve beyond e max , the po- 
tential lepton number would be completely depleted when e 
reaches e c .o. ~ 0.987 (i.e., £(e c .o.) = 0). 



other neutrinos have rates which are slow compared to 
the expansion rate of the universe. This means that these 
processes will be unable to redistribute effectively the 
active neutrino occupation numbers and so they cannot 
morph the v a distribution into a thermal distribution. 
This has consequences for the lepton capture rates on 
nucleons as we will discuss below in section IV. 

In either the (unphysical) smooth and continuous res- 
onance sweep scenario or in some stochastic resonance 
sweep case both the active neutrino and resulting ster- 
ile neutrino distribution functions will be non-thermal in 
character. Since we cannot solve for resonance sweep be- 
yond e ma x we do not know the final active and sterile 
neutrino energy spectra. However, for the purposes of 
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FIG. 3: The values of e c . . (solid line) and e max (dashed line) 
are shown as functions of total initial potential lepton number 
in the limit of a smooth and continuous resonance sweep and 
with the assumption that full active neutrino equilibration 
obtains (L Vc = L Vfi = L„ T ). 



constraints and general guidelines, we will find that the 
idealized smooth and continuous resonance sweep sce- 
nario provides the basis for lower limits on the effects of 
neutrino spectral distortion. 

It is useful to consider the solution for the cut-off en- 
ergy e c . . and the peak energy e max in the smooth and 
continuous resonance sweep case. To get the first of these 
quantities, we force the system to evolve continuously 
(e.g., along the solid curve in Fig. [21 all the way to com- 
plete lepton number depletion and solve 



£ 



initial 



l 



4 F 2 (0) 



1 



-dx. 



(27) 



or £(e c . .) = 0. The second of these quantities is the 
solution of Eq. (|26p. Both of these solutions are shown 
as functions of initial potential lepton number, £ lnltlal in 
Fig. (J3J . In this figure it is assumed that the active neu- 
trinos are fully "equilibrated" initially (before any flavor 
transformation) with L Ve = L v = L Ur . 

It is obvious from Eq. 1)27(1 that there is a maximum 
value of the initial potential lepton number for which a 
solution is obtainable when v a ^ v s is the only operative 
neutrino flavor conversion channel. This maximum is 
given by the limit where e c ., 



oo, 



C 



initial 



3 F 2 (i lva ) 

4 F 2 (0) ' 



(28) 



Scenarios where the bulk of the initial potential lepton 
number is contained in seas of another flavor of active 
neutrinos may not allow v a ^ v s conversion to leave 
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a zero final potential lepton number. This is a simple 
consequence of the post Weak Decoupling conservation of 
numbers of neutrinos of all kinds. Of course, active-active 
neutrino flavor transformation in the channels v a ^ vp 
(a, (3 = e, (a, t) can alter this picture significantly and will 
be discussed below. 



C. Efficiency of Neutrino Flavor Conversion: 
Adiabaticity 



proper, physical entropy density in radiation-dominated 
conditions is S ~ (2ir 2 / 45) g S T 3 , where g s is closely re- 
lated to g and we can take g s s» g. We can take the co- 
moving volume element to be the cube of the scale factor 
a in the Friedman-Lemaitre-Robertson- Walker (FLRW) 
metric, so that a 3 S is invariant with FLRW time coor- 
dinate t and therefore g 1 / 3 aT is constant. In turn, this 
implies that the fractional rate of change of the temper- 
ature is related to the expansion rate and the fractional 
rate of change of the statistical weight by 



From the previous discussion it is clear that the effi- 
ciency of active-sterile neutrino flavor conversion at MSW 
resonances is a key issue in resolving how these reso- 
nances sweep with scaled energy e. This is especially 
true for values of scaled resonance energy beyond e max , 
where we argued that if the resonance condition can be 
met neutrino flavor transformation was likely to be ef- 
ficient for the typical neutrino mass/mixing parameters 
we consider here. It is the adiabaticity of neutrino prop- 
agation which determines transformation efficiency both 
in the active-sterile and active-active channels. We point 
out here that adiabaticity parameters are high for our 
chosen epoch and neutrino mass/mixing characteristics, 
essentially because these parameters are proportional to 
the ratio of a gravitational time scale (the causal horizon) 
to a weak timescale (in-medium oscillation time). 

The causal horizon (particle horizon) is the proper dis- 
tance traversed by a null signal in the age of the universe 
t. In radiation-dominated conditions in the early universe 
this is (setting c = 1) 



dn (f) = 2t = H- 1 , 
where the local Hubble expansion rate is 



H 



37T 
~90" 



1/2 



,1/2 



TOpl 



(29) 



(30) 



Here m pl » 1.221 x 10 22 MeV is the Planck mass. The 
statistical weight for a relativistic boson species i is (fffeL, 
while that for a relativistic fermion species j is (gf)j- 
These are related to the total statistical weight g by a 
sum over all particle species i and j with relativistic kine- 
matics with equilibrium or near equilibrium energy dis- 
tribution functions and energy densities in the plasma at 
temperature T, given by g = £\ (g b ) i + (7/8) £V (g f ) r 
In the epoch between Weak Decoupling and Weak Freeze 
Out and BBN, photons, e^- pairs and the active neutri- 
nos are relativistic and appreciably populated so that 
g w 10.75 and t m (0.74 s) (10.75/.g) 1/2 (MeV/T) 2 . The 
spectral distortions and extra energy density stemming 
from the net lepton numbers considered in this paper are 
usually small effects, causing deviations of the expansion 
rate from that given above by less than a few percent in 
many cases. 

Homogeneity and isotropy in the early universe imply 
that the entropy in a co-moving volume is conserved. The 
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3H 



(31) 



At lower temperatures, where the thermal potential 
can be neglected, the potential governing neutrino fla- 
vor transformation is the difference of the Hamiltonians 
(e.g., Eq.s Q , JSJ , © , 10 ) for the transforming neutrino 
species. For the active-sterile channel v a ^ i/ s} for ex- 
ample, we have V » H(v a ) — H(v s ). The appropriate 
potentials for the active-active neutrino flavor transfor- 
mation channels follow in like manner. 

The density scale height for the early universe depends 
on the neutrino flavor transformation channel and is de- 
fined as 



H 



1 dV 
V~dt 



(32) 



1 



1 



9/9 
3H 



L/C 



3H 



Here the approximation on the second line is for active- 
sterile neutrino flavor transformation and follows on ne- 
glecting the thermal potential C and using Eq. 131fl . 
When the statistical weight and the lepton numbers are 
not changing rapidly the density scale height is roughly a 
third of the horizon scale. This is ~ 10 km at the epoch 
we are considering here. 

Define A = 6m 2 /2E V . It can be shown that the ratio 
of the difference of the squares of the effective masses in 
matter to twice the neutrino energy is 



A, 



Sm 2 



off 



2E V 



(33) 



(A cos 26 - V) 2 + (A sin 26 + B eT ) 2 , 



where 6 is the appropriate effective two-by-two vacuum 
mixing angle and where V = A + B + C is the appropri- 
ate potential for the transformation channel. Here B eT 
is the flavor-off diagonal potential as defined in Qian & 
Fuller 1995 |2(|. The flavor basis off-diagonal potential 
vanishes, B eT = 0, for any active-sterile mixing channel. 

The effective matter (in-medium) mixing angle 6m for 
a neutrino transformation channel with potential V and 
effective vacuum mixing angle 6 satisfies 
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sin 2 29m = 



A 2 sin 2 20 (l + 2E u B eT /Sm 2 sin 29) 2 
(A cos 26> - V) 2 + A 2 sin 2 20(1 + 2E u B eT /Sm 2 sin 20) 2 ' 



(34) 



r 



The effective matter mixing angie for the antineutrinos 
in this channel, 9m, satisfies a an expression which has 
opposite signs for the potentials B, A, and B eT , but which 
is otherwise identical. 

The change in the potential required to drop the ef- 
fective matter mixing from the maximal resonant value 
(9m = 7r /4) to a value where sin 2 29m = 1/2 is termed 
the resonance width and is 



SV ~ A sin 29 



1 



Sm 2 sin 29 



(35) 



The physical width in space, or in FLRW coordinate time 
t, corresponding to this potential width is 



Htan20 



1 dV 


^SV 


V~dT 





(36) 



1 



2E u B e 



Sm 2 sin 29 

The local neutrino oscillation length at resonance is 
4tvE v 2ir 2ir 



t res 

osc 



6 mi 



icir 



SV 



(37) 



where the latter approximation is good only at resonance. 
We can define the dimcnsionless adiabaticity parameter 



as proportional to the ratio of the resonance width and 
the neutrino oscillation length at resonance: 



2tt 



St 

J res 

osc 



StSV 



(38) 



Sm 2 H sin 2 29 



2E V 



cos 20 



2E V B 



Sm 2 sin 29 



This parameter can be evaluated anywhere in the evolu- 
tion of neutrino flavors, even well away from resonances 
and it will serve to gauge the degree to which neutrinos 
tend to remain in mass eigenstates. The Landau-Zener 
jump probability, assuming a linear change in potential 
across the resonance width, is Plz ~ ex p(~ 7r 7/2), so 
that it is clear that a large value of the adiabaticity pa- 
rameter corresponds to a small probability of jumping 
between mass eigenstate tracks and, hence, efficient fla- 
vor conversion at asymptotically large distance (many 
resonance widths) from resonance. 

Folding in the expansion rate in radiation-dominated 
conditions, using the conservation of co- moving entropy 
density, and assuming that we can neglect the thermal 
potential C, we can show that the adiabaticity param- 
eter for neutrino propagation through an active-sterile 
resonance is 



J 



V5C 3/4 (3) (W) 1/4 m pl G| /4 



2 1 /8 7r 3 
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' Sm 2 ' 




.leV 2 . 



1/2 
1/4 
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r £ 3/4- 




' sin 2 2(9 " 


[ e l/4 




cos 7 /4 29 


" C ' 


3/4 




9/9 
3H 


0.01 







c/c 



g/g 
m 
-i 



L/C 



3H 



3H 



sin 2 26» 
1.77 x 10" 



r 



(39) 



In these expressions e = E v jT is the scaled energy of a 
neutrino at resonance in a channel v a ^ v s character- 
ized by the difference of the squares of the appropriate 
vacuum mass eigenvalues, Sm 2 . It is obvious from these 
considerations that neutrino flavor transformation will be 
efficient at resonance (i.e., 7 3> 1) over broad ranges of 
energy for the regime of the early universe between Weak 
Decoupling and Weak Freeze Out even for very small ef- 
fective vacuum mixing angle 9. 

Eq. Ij39(l shows that two trends can eventually destroy 
adiabaticity and, therefore, large scale resonant active- 
sterile neutrino flavor transformation. As active neutri- 



nos are converted C is reduced and this reduces 7. In 
turn, the fractional rate of destruction of C compared 
with the Hubble parameter can be become significant, 
especially if C is small, and this can also reduce 7. 



D. Active-active Neutrino Flavor Conversion and 
Equilibration 

Active neutrinos (v e ,v e ,v tl ,v ill v ri v T ) transforming 
among themselves on time scales comparable to or 
shorter than that of the active-sterile conversion chan- 



10 




> 



artificial 
1 eV 2 offset 

I 



e=i 

L Ve =0.0343 



1.20 



T (MeV) 



FIG. 4: Level crossing diagram for the case with lepton num- 
bers as shown and for scaled neutrino energy e = 1. The vac- 
uum mass-squared eigenvalue for the (mostly) sterile state is 
taken as m 2 = 10 eV 2 . This is shown as the dashed curve 
labeled v s . An artificial (exaggerated) 1 eV 2 offset between 
the vacuum mass-squared eigenvalues m 2 . and m% has been 
added so that the curves labeled with v* and v* are sepa- 
rated for clarity. In reality, the top curve should be split from 
the lower curve by 6m 2 ~ 3 x 10~ 3 eV 2 . Conversion in the 
channel i/ e — > v s is as described in the text. 



nel can alter significantly the scenario for sterile neu- 
trino production given above. This is apt to be the case 
if active-active neutrino mixing in medium is large and 
efficient over a broad range of neutrino energies. Active- 
sterile neutrino flavor conversion tends to be slow because 
it occurs through MSW resonances and the rate at which 
these resonances sweep through the neutrino distribution 
functions is determined by the expansion of the universe, 
a slow gravitational time scale. 

Coherent neutrino flavor conversion in active-active 
channels in the early universe can be dominated by the 
flavor off-diagonal potential. Large in-medium mixing 
angles can accompany the synchronization seen in cal- 
culations of active-active mixing in supernovae and the 
early universe [l7| . If active-active neutrino flavor trans- 
formation is efficient, then lepton numbers in different 
active neutrino species can be quickly equilibrated, mean- 
ing instantaneous equal lepton numbers. 

The flavor diagonal neutrino forward scattering poten- 
tial in an active-active channel v a vp is A + B = 
H (v a ) — H (vp). If there is an initial disparity in lep- 
ton number in these two flavors then matter-enhanced 
or -suppressed transformation will go in the direction of 
reducing this disparity. Though initially the flavor off- 



diagonal potential B eT w 0, as soon as flavor transforma- 
tion begins this potential comes up. 

The interplay of matter-enhanced coupled active- 
sterile and active-active neutrino flavor transformation 
can be complicated and difficult to follow numerically. 
The size of the debit in the v e or v e distributions {e.g., 
the final value of e c . . or e max in a continuous sweep sce- 
nario) may be much more complicated in the general 4x4 
case than the scenario outlined above for "simple" 2x2 
Va v s interconversion. We can, however, identify a few 
cases where we can at least outline the course of neutrino 
flavor conversion as the universe expands and cools. We 
will therefore consider two limits: (1) no active-active 
mixing; and (2) efficient active-active mixing that guar- 
antees that lepton numbers in active species are always 
the same (instantaneous equilibration). 



E. Inefficient Active- Active Neutrino Flavor 
Conversion 

Consider first the case where we neglect active-active 
neutrino mixing effects. In this case we could have 
initial lepton numbers that are not fully equilibrated. 
For example, we could have a scenario where initially 
L Vii < L u — L Vt . In this case the and v T experi- 
ence the largest effective potentials, and hence have the 
largest effective masses at a given temperature (epoch) in 
the early universe. Therefore, the first (highest temper- 
ature) resonance occurs for v s with z/ e , as illustrated in 
Fig. This resonance will destroy lepton number, as 
will the subsequent u* # v s resonance, and it will leave 
a distorted v e spectrum. 

Here we follow Ref. 27] and define linear combinations 
of the muon and tauon neutrino flavor states 



K) 



V2 

~ V2 ' 



(40) 
(41) 



This reduces the 4x4 mixing problem of three active neu- 
trinos and a sterile neutrino into a 3 x 3 problem with \ v*) 
decoupled (a mass eigenstate in vacuum with no mixing 
with the other neutrinos). This reduction in dimensional- 
ity of the neutrino mixing problem works in vacuum only 
if the muon and tauon neutrinos are maximally mixed. 
It will be valid in medium only if, additionally, these two 
neutrino flavors experience identical matter interactions. 
This latter condition is met if L v i = L Vt . This sym- 
metry condition will be respected so long as muon and 
tauon neutrinos behave and transform identically. In- 
deed, the second resonance encountered as the universe 
cools, is* v s , respects this condition as the \u*) state 
consists of equal parts muon and tauon states. 

The sterile neutrinos produced through the v e — > v s 
resonance are subsequently transformed into v* at the 
second resonance at lower temperature, as depicted for a 
particular set of initial lepton numbers in Fig. This 
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resonance also converts the v* into the sterile state, so 
that the final abundance of sterile neutrinos results from 
the conversion of neutrinos which were originally in the 
Ufj, and v r distributions. Since these distributions have 
higher lepton number than resides in the v e /ve seas, the 
final number density of sterile neutrinos will be larger 
than the number of v e missing from the ^-distribution. 
As we will see, this case may be more likely to be in 
conflict with massive neutrino dark matter constraints. 

If we temporarily ignore the effect of active-active neu- 
trino flavor transformations, then we can make some gen- 
eral statements about the change in the lepton numbers 
and e for the v e or v e distributions in this case of unequal 
L„ c and L v — L Ut . If the potential lepton number for 
electron flavor neutrinos C e is driven to zero first then the 
changes in the individual active neutrino lepton numbers 
must be related by 

2AL, s + AL„» = -4 nitiaJ , (42) 

where here L u * = [L Vll + L Ur ) /2, and AL^ = An„ e /n 7 
and AL Ull = An Ufi /n 7 . 

Of course, if C e is driven to zero before £* , then con- 
version of ^ e 's at the first resonance will cease while con- 
version in the channel i>* — > v s continues until £ M » is 
reduced to zero. This will leave C e < which will result 
in anti-electron neutrino transformation v e — > v s , leaving 
a non-thermal deficit in the v & distribution. 

This is likely temporary, however. The v e potential 
is zero at the point where C e first vanishes. Thereafter, 
with reduction in , the v e potential's magnitude first 
increases, but then decreases as the universe expands and 
the temperature drops. This can be seen from Eq. ITU 
and on noting that the potential behaves like V ~ C e T 3 . 
The evolution of the potential with time is determined 
by the competition between two effects. The conversion 
of lepton number discussed above makes C e more neg- 
ative and larger in magnitude, while the expansion of 
the universe decreases T. Therefore, the antineutrinos 
could experience two resonances: (1) first when C e be- 
comes sufficiently negative that the potential V becomes 
large enough in magnitude to satisfy Eq. lfH|) ; and (2) 
subsequently when the temperature drops enough that 
this condition is again satisfied. At the first resonance 
we have v e — > D Sl but these steriles are reconverted at 
the second resonance, v s — > v e . This can also be viewed 
from the perspective of resonance sweep. Note that e res 
is infinite when C e first crosses zero, but then decreases 
as C e becomes larger in magnitude as lepton number is 
converted, but then "turns around" and begins to sweep 
toward higher energy again once T becomes low enough. 
This process is directly analogous to the re-conversion 
of v s neutrinos in neutrino-heated outflow in supernovae 

H 

In this scenario it is the mu and tau neutrinos, ulti- 
mately, that are converted to sterile neutrinos so that 
the numbers and kinds of converted active neutrinos are 



given by 

~ A ( A n,. + An'„. ) - — (An„ e + An Pe ) , (43) 
^ n 7 V M *v n 7 

where An,,* and An', are the number of v* neutrinos 
converted before and after C e first vanishes, respectively. 
Likewise, An„ e electron neutrinos are converted before 
C e first vanishes and An Pc electron antineutrinos after- 
ward, though these P e 's are eventually returned to the 
distribution. 

There is an additional complication: in the case that all 
three lepton numbers are equal, the mass-squared differ- 
ences between the active states are approximately given 
by their vacuum values, which are quite small. The three 
resonances depicted in Fig. will then be very close to- 
gether, and in fact may overlap if the resonance width is 
sizable. 

If the resonances do not overlap, the lepton number 
destroying resonance will take place between v\ and v Sl 
where v\ is the lightest neutrino mass eigenstate. Since 
v\ has a large v e component, this will leave a non-thermal 
v e distribution, and in addition there will be smaller non- 
thermal distortions of the and v r spectra. In the case 
that the resonances do overlap, the full details of the evo- 
lution will be quite complicated, but a similar outcome 
is obtained nonetheless. To summarize, in all cases a non 
thermal v e spectrum results. 

F. Efficient Active- Active Mixing: Instantaneously 
Equilibrated Lepton Numbers 

Let us now consider the limit where in addi- 
tion to the active-sterile MSW transitions, oscilla- 
tions/transformations between/among the three active 
neutrinos occur simultaneously and are efficient. If 
active-active mixing among all the active flavors is instan- 
taneous and efficient then we only need to consider the 
case where the lepton numbers are equal, L e = = L T , 
both initially and as active-sterile transformation pro- 
ceeds. It has been shown that large angle mixing between 
the three active neutrino species results in the system 
being driven toward such an equilibrated state |l7j at a 
temperature of T > 2 MeV. 

An obvious additional effect of efficient active-active 
oscillations will be to partially refill any hole that was left 
in the v e distribution. It is important to note, though, 
that this refilling cannot be complete. For maximal 
Ve. — JVt mixing, the hole in the distribution can be only 
partially refilled. In vacuum the measured solar neutrino 
mixing angle is less than maximal, d SQ \sa — 32.5°, and 
U e 3 is relatively small. In medium, at best we will ob- 
tain maximal matter mixing angles in the limit where the 
flavor off-diagonal potential is large. We again therefore 
expect about to 2/3 refilling at most, so that a non- 
thermal v e spectrum is always obtained by the epoch of 
Weak Freeze Out. Even if it were somehow possible for 
the resonance to effectively involve only v s and v^/v-r, 
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FIG. 5: Final active neutrino energy distribution function 
f Ua for Cases 1 (short dashed line), 2 (dot dashed line), and 
3 (solid line) for the particular case of a continuous resonance 
sweep scenario and in the instantaneous active-active mixing 
limit as described in the text. Here a = e,^i,r: all species 
have the same distribution function. The long dashed line 
and its continuation as a solid line shows the original thermal 
distribution function common to all active flavors. The par- 
ticular scenario shown here has L Vs = L Vt = L u = 0.1, so 
that ei « 3.8, £2 ~ 2.3, and £3 « 1.85 



active-active oscillations would again act to refill the hole 
in the resulting non thermal and/or v T spectra, and 
in so doing create a non-thermal v e distribution. 
We can identify three cases. 

Case 1: We have only one sterile neutrino species and 
only one channel for its production, v a — > v s . If this 
channel is, e.g., v e — > v Sl then the neutrinos in the v T and 
distributions will, in the limit of instantaneous max- 
imal mixing, partially fill in the hole left by the active- 
sterile conversion process. Given the boundary condition 
of equal lepton numbers in all active flavors at all times, 
a continuous smooth resonance sweep scenario will leave 
each active neutrino distribution with a low energy "hole" 
with 2/3 of the normal population out to some value of 
scaled neutrino energy t\ . In terms of the initial potential 
lepton number £ mlt = C e = £ M = C T , this is obtained 
by solving the integral equation 



2C 



init 



1 



F 2 (0) 



x dx 

x-ri _|_ I ' 



(44) 



Of course, the continuous resonance sweep scenario is 
an unphysical idealization and the actual energy spec- 
tra will likely be far more complicated, as argued above. 
However, no matter the resonance sweep physics, in this 
instantaneous mixing limit the numbers (number densi- 
ties) of active neutrinos in each flavor are equal and their 
energy spectra are identical. In the continuous resonance 
sweep scenario, in the scaled energy interval to e\, the 
deficit of neutrinos, An„, is 1 — 2/3 = 1/3 of the original 
population and this deficit is identical for each active fla- 
vor. The entire original population of one of the active 
neutrino flavors in this scaled energy interval is converted 
to sterile neutrinos so the number density of steriles will 
be n s = 3An„ and this is related to the initial potential 



3An„ _ ^£init 

n 7 4 



(45) 



If there is one light sterile neutrino, there may be oth- 
ers. In fact it has been claimed that two sterile species 
are a better fit to the LSND data than just one [2£| ■ So 
this suggests Cases 2 and 3. 

Case 2: Allow two channels of sterile neutrino produc- 
tion and two kinds of light sterile neutrinos v s \ and v S 2- 
As an example, we could have v e — > v s \ and — > v S 2, 
but again with instantaneous mixing among all the active 
neutrino flavors. A continuous resonance sweep scenario 
will leave each active neutrino distribution with a low en- 
ergy hole now with 1/3 of the normal population out to 
some value of scaled neutrino energy £2. In terms of the 



initial potential lepton number C 



init 



C = C 



this is obtained by solving the integral equation 
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init 
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F 2 (0) 



e x-r, + 1 
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(46) 



Again, the numbers of active neutrinos in each flavor in 
scaled energy interval to £2 are equal and so are the 
deficits, An u , which are now 2/3 of the original popula- 
tion in this interval. Now, however, the original popula- 
tions of two active flavors in this interval are converted 
to sterile species so the total number density of sterile 
neutrinos of all kinds is n s = 2 • (3/2) An^ = 3Ani, and 
we have the same relation between total sterile neutrino 
number density n s , deficit per flavor An„, and initial po- 
tential lepton number as in Eq. (|45() . 

For a given initial potential lepton number, this is the 
same total number of sterile neutrinos (of all kinds) pro- 
duced as in Case 1. However, since there are now two 
channels for v s production, e 2 is smaller than ei. In Case 
1, ei is relatively larger because as ^ e 's are converted to 
steriles two active neutrino distributions compensate by 
feeding neutrinos into the hole, forcing the resonance to 
sweep further (higher in energy) through the v e distri- 
bution to erase the net lepton numbers. In Case 2 only 
one active neutrino distribution remains to compensate 
for the hole. 

Case 3: Allow all three active neutrinos to convert 
simultaneously to three kinds of light sterile neutrinos 
v s i, Vs2-, ^s3- Now a smooth resonance sweep scenario 
will leave each active neutrino distribution with a low 
energy hole with zero population out to some value of 
scaled neutrino energy £3. In terms of the initial potential 
""' — r - Cn = C T , this is obtained 



lepton number £ = C e = 
by solving the integral equation 



^ /-init * 

3 ~ F 2 (0) Jo 



e x-i) + I 



(47) 



Now the deficits An„ in each active flavor are equal to the 
original active neutrino populations in the scaled energy 
interval to £3. Since all three active species transform 
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to sterile neutrinos, the total number density of steriles 
of all kinds is n s — 3An„. The relation between n s , 
An„ and the initial potential lepton number is the same 
as in Eq. (|45|l . For a given £ lmt this is the same total 
number of sterile neutrinos produced as in Cases 1 and 
2. However, our three cases will have ei > 62 > £3 for 
a given initial potential lepton number, for the reasons 
indicated in the last paragraph. 

For Cases 1, 2, and 3 the active neutrino distribution 
functions will be left with population deficits relative to 
the thermal case. This is shown in Fig. (JSJ for the partic- 
ular scenario where each active flavor starts out with lep- 
ton number L Ve = L Ut = L u =0.1. Solving the above 
equations for the three cases yields t\ « 3.8, £2 ~ 2.3, 
and 63 ss 1.85 in this example. 

In obvious fashion all of the above discussion applies to 
v s production if the initial lepton numbers are negative. 
We should also note that the actual active and sterile 
neutrino energy distributions in all of the limits consid- 
ered here may differ considerably from those shown in the 
figures. This is partly because the continuous resonance 
sweep cannot continue to completion as described in the 
last section but may skip to higher scaled energy dis- 
continuously. However, another source of difference from 
the simple spectra shown in the figures may be because 
multiple neutrino mixing can be complicated. 

We have here presented a picture where the sterile neu- 
trino undergoes a resonance with one of the active neu- 
trino flavors, v a . However, if the initial lepton numbers 
are equal, the resonance will instead occur between v 8 
and a superposition of the three active neutrinos. In 
principle, all three of the active neutrinos may mix with 
the sterile, so the MSW resonance which is responsible 
for lepton number destruction may occur for the sterile 
neutrino and a linear superposition of the three active 
neutrinos. For example, in the so called "3+1" LSND in- 
spired mixing scheme, both #14 and #24 are required to be 
non-zero, and the sterile effectively mixes with all three 
active neutrinos (as there will also be indirect v T — v a 
mixing) . 

In any case, post-decoupling neutrino mixing cannot 
completely undo spectral distortions. We conclude that 
a sterile neutrino in the mass range of interest is almost 
certain to leave non-thermal active neutrino distribution 
functions if the lepton number is significant. 



III. CONSTRAINTS ON STERILE NEUTRINOS 
AND LEPTON NUMBERS 



The entire plausible range of sterile neutrino masses 
and net lepton numbers of interest is not likely to be 
consistent with all of the current observational bounds. 
For example, we may demand that the initial net lepton 
numbers are large enough to suppress the production of 
fully thermalized seas of v s and v s . Eq. @ shows that the 
lepton number necessary for suppression of thermal ster- 
ile neutrino production depends both on neutrino mass 



and neutrino average energy. We will hold off on consid- 
ering BBN effects/limits until the next section. 

As discussed in the introduction, a population of sterile 
neutrinos could provide a sufficiently large contribution 
to the dark matter density, depending again on sterile 
neutrino mass, to run afoul of large scale structure/CMB 
bounds 8J. The above-cited analysis of the SDSS data 
|lOj using CMB anisotropy limits, galaxy clustering and 
bias, and coupled with the matter power spectrum in- 
ferred from the Lyman-alpha forest suggest a limit on the 
neutrino mass of 0.79 eV (95% CL). This corresponds to 
a limit on the neutrino closure fraction 



n^h 2 < 0.0084 (95% CL), 



(48) 



where h is the Hubble parameter at the current epoch 
in units of 100 km s" 1 Mpc~ x . This is comparable to the 
older Wilkinson Microwave Anisotropy Probe (WMAP) 
bound, fl^h 2 < 0.0076 (95% CL). However, the 
Eq. (|48|l bound is more appropriate here as it assumes 
a "3+1" neutrino mass scenario in contrast to the three 
neutrinos with a common mass assumed in the WMAP 
analysis. Of course, this mass limit and our inferred limit 
on the closure fraction are approximate because our ster- 
ile neutrinos have nonthermal energy spectra. Adopting 
the Eq. I)48|) bound suggests that thermal distributions 
of v a and v a neutrinos are acceptable only if they have 
rest masses 



< 0.79 eV 



2F 2 (0) 



F 2 {r, Va )+F 2 (-r) Va ) 



0.0084 



, (49) 



where a = e, /1, r, s. We can connect this with the 
a putative thermal sterile neutrino sea by noting that 

m„ a ~ (Sm 2 ^ 1 ^ 2 . So, for example, <5m 2 s > 0.63 eV 2 is 
disallowed if all the sterile neutrino species have thermal 
distributions. This would eliminate much of the LSND- 
inspired sterile neutrino mass range. 

However, the coherent sterile neutrino production sce- 
narios discussed above may do better at creeping in un- 
der the closure contribution bound. For one thing, only 
v s (or v s ) and not its opposite helicity partner are pro- 
duced coherently. Furthermore, the sterile neutrinos are 
produced in numbers of order the initial lepton number. 
This will be smaller than a general thermal population. 

At the epoch of coherent sterile neutrino production 
the ratio of the number of active neutrinos An^ con- 
verted to steriles to the number total density of a ther- 
mal distribution of v a plus v a neutrinos is in the ratio of 
the closure contributions of a sterile species to thermal 
neutrino species: 



n,h 2 



R, 



N„Ar 



(50) 



where N s is the number of active-sterile mixing channels 
operating in the production of sterile neutrinos. In turn 
it can be shown that in the continuous resonance sweep 
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FIG. 6: Constraints on the the ranges of active-sterile mass- 
squared splitting and potential lepton number as derived in 
the smooth and continuous resonance sweep limit. Parame- 
ter ranges which give sufficient sterile neutrino production to 
exceed the bound on the neutrino closure fraction are shown 
cross hatched, as are parameter ranges which allow for com- 
plete or nearly complete thermal, undiluted energy distribu- 
tion functions for a sterile species. The upper solid line is for 
efficient active-active mixing in Cases 1, 2 , or 3, while the 
upper dashed line gives the constraint for v a — > v s with no 
active-active mixing. 



scenario 



l 



F 2 {0) 



*2(0) 



(51) 

where e c . . and the degeneracy parameter r\ are values 
consistent with the particular sterile neutrino production 
scheme. From these relations we can show that 



n s h 2 fa (1.062 x 10" 



(C 










[eV 2 \ 



1/2 



(52) 



where C is an appropriate potential lepton number and 
where (3 is a parameter that is related to the particu- 
lar sterile neutrino production scheme and the number 
of active-sterile channels in that scheme. For example, 
j3 = 2 for Cases 1, 2, and 3 of the efficient active-active 
limit, whereas (3 = 4/3 for v a — > v s only with no active- 
active mixing. All of these constraints are summarized 
in Fig. ©. 

Though we have employed the continuous resonance 
sweep scenario, the constraints shown stem from closure 
fraction considerations and, therefore, depend principally 
on numbers of sterile neutrinos produced and on their 



rest mass. The number of sterile neutrinos is tied to the 
initial potential lepton number and is insensitive to the 
details of resonance sweep physics when the all the lep- 
ton number is destroyed and the number of neutrinos of 
all kinds is fixed. We expect the latter condition to be 
a good approximation post Weak Decoupling. The for- 
mer assumption is a good one over most of the range 
of neutrino masses/mixing considered here because neu- 
trino flavor evolution is highly adiabatic for these param- 
eters as shown above. 

However, the actual energy spectra of the sterile neu- 
trinos may come into play as precision in observations of 
the matter fluctuation power spectrum increase. Current 
constraints are most strongly dependent on the suppres- 
sion of small scale power. Future constraints may be able 
to constrain the collisionless damping scale of neutrinos — 
at large scales — much more strongly, and are most im- 
portant in regard to spectral constraints. The collision- 
less damping scale for neutrinos is essentially their free 
streaming length, which of course depends on the neu- 
trino velocities which, in turn, depend on the neutrino 
energy spectra. 



IV. NON-THERMAL NEUTRINO ENERGY 
SPECTRA: ALTERATION OF THE 
RELATIONSHIP BETWEEN LEPTON NUMBERS 
AND PRIMORDIAL 4 HE ABUNDANCE 



Primordial nucleosynthesis is a freeze out from Nu- 
clear Statistical Equilibrium (NSE). In NSE the abun- 
dance of nuclei is set by a competition between disorder 
(entropy) and binding energy. In the early universe al- 
pha particles win this competition because the entropy 
per baryon is very high and alpha particles have a bind- 
ing per nucleon not terribly different from iron. At 
a temperature T a ~ 100 keV, alpha particles form ag- 
gressively, incorporating essentially all neutrons (all but 
~ 1/10 5 ). Therefore, the primordial 4 He yield is de- 
termined roughly by the neutron-to-proton ratio n/p at 
T a . In mass fraction, this is X a ps 2(n/p)/(n/p + 1), or 
25% for n/p = 1/7. The standard BBN 4 He mass frac- 
tion yield prediction is (24.85 ± 0.05)% using the CMB 
(Cosmic Microwave Background) anisotropy-determincd 
baryon density @. (The baryon closure fraction as de- 
rived from the deuterium abundance [3fl | is consistent 
with the CMB-derived value.) 

The observationally inferred primordial helium abun- 
dance has a long and troubled history. One group pegs 
this abundance at 0.238±0.002±0.005 HJ, while another 
using similar but not identical compact blue galaxy data 
estimates 0.2421 ± 0.0021 These values are quite 

restrictive. However, these older estimates may now be 
superseded by more recent analyses as discussed in the 
Introduction. In particular, a more detailed analysis of 
the helium and hydrogen emission lines done in Ref. Q 
suggests that the allowable range of mass fraction for 
primordial 4 He is 0.232 to 0.258. This is fairly generous 
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compared to previous "limits," but it is still a good bet 
that a 5% or 10% increase in the calculated, predicted 
yield in 4 He would be an unwelcome development. 

The relationship between neutrino physics and/or lep- 
ton numbers and the primordial helium abundance re- 
mains a cornerstone of modern cosmology. Distortions 
in the v e and v e energy spectra stemming from active- 
sterile conversion can alter lepton capture rates on nu- 
cleons and thereby change n/p and the 4 He yield com- 
pared to a standard BBN scenario with thermal neutrino 
energy spectra. 

The neutron-to-proton ratio is set by the competition 
of the expansion rate of the universe and the rates of the 
following lepton capture/decay processes: 

u e + n # p + e~, (53) 



Ve + P n + e + , (54) 



(55) 



We denote the forward and reverse rates of the first pro- 
cess as \ VcTl and A e - p , respectively. Likewise, the forward 
and reverse rates of the second process are A PcP and A e +„, 
respectively, while those of the third process are A„ decay 
and A pePe , respectively. 

The lepton capture processes' influence on the isospin 
state of nucleons can be appreciable, even for the post 
Weak Decoupling epoch. This is because the number 
densities of relativistic neutrinos and charged leptons are 
some 10 orders of magnitude larger than the baryon den- 
sity. At high enough temperature (T 3> lMeV), where 
these rates are very fast, the isospin of any nucleon will 
flip from neutron to proton and back at a rate which is 
rapid compared to the expansion rate of the universe, 
establishing a steady state equilibrium. 

As the universe expands and the temperature drops 
the rates of the lepton capture processes will drop off 
quickly. Eventually the lepton capture rates will fall be- 
low the expansion rate and n/p will be frozen in, save for 
free neutron decay. Traditionally, this "weak freeze-out" 
epoch is taken to be T w f D » 0.7 MeV. 

However, there is no sharp freeze-out of isospin. In 
fact, the neutron-to-proton ratio n/p is modified by the 
lepton capture reactions down to temperatures of several 
hundred keV and by neutron decay through the epoch of 
alpha particle formation T a . 

The evolution of the electron fraction Y e = 1/(1 + n/p) 
throughout the expansion is governed by 



dY e 
dt 



A, 



Y K At 



(56) 



where the sum of the rates of the neutron destroying 
processes is A„ = \ UaTL + A e +„ + A Ilc i C cay and the sum of 
all weak isospin changing rates is A tot = A n + A PcP + 
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FIG. 7: The neutron to proton ratio as a function of tem- 
perature for the indicated lepton numbers. The lower solid 
curve gives the case with lepton numbers alone with no active- 
sterile neutrino conversion. The dashed line is the equilib- 
rium prediction for the n/p ratio in this case. The upper 
solid curve is the n/p ratio for the same lepton numbers, 
but now with active-sterile conversion in the v e — * v s chan- 
nel with <5m 2 = 10 eV 2 , while the middle solid curve is for 
Sm 2 = leV 2 . The upper dotted curve is the steady state 
equilibrium prediction for the n/p ratio for the Sm 2 = 10 eV 2 
case, while the lower dotted curve gives the same prediction 
1 eV 2 case. 



In the limit where the isospin flip rate is fast compared 
to the expansion rate H , the neutron-to-proton ratio has 
a steady state equilibrium value given by |3Jfl 



n 
P 



Ay eP + A, 



-p ~t~ Ap e iy e 



Aiy e n ~t~ A e + n -|- An decay 
Ap e p ~\~ A e — p 



(57) 



Aiy e n ~t~ A e + n 



This solution corresponds to the fixed point dY e /dt = 
in Eq. I|56() . where Y e — A„/A to t- The second approxi- 
mation in Eq. (|57|l is valid at temperatures high enough 
that the rates of the three body processes can be ne- 
glected relative to the lepton capture rates. 

As the universe expands and the temperature drops, 
the relative values of these rates change and, hence, so 
does the neutron-to-proton ratio. Certainly at tempera- 
tures T > T w f the three body lepton capture and free 
neutron decay processes have rates which are unimpor- 
tant compared to those of the lepton capture rates. For 
temperatures T ^> T w f Q , typical lepton energies are large 
compared to the energy thresholds in the forward rate 
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of the process in Eq. 1(54(1 and the reverse rate for the 
process in Eq. (|53|l . so that if the lepton numbers are 
small we would have n/p ~ 1, or Y e w 1/2. During 
the epoch T w f Q > T > T a the lepton capture processes 
gradually give way to free neutron decay as the princi- 
pal n/p-altering mechanism. Over this time period e^- 
annihilation heats the photon/electron plasma relative to 
the neutrinos, further altering the weak rates (including 
the free neutron decay) by modifying the neutrino and 
antineutrino distribution functions relative to those for 
e ± . (See Ref. 0| for a discussion of primordial nucle- 
osynthesis.) 

If the electron neutrinos and antincutrinos and the 
electrons and positrons all have Fermi-Dirac energy spec- 
tra, then Eq. (|BTjl can be reduced to [sE 

n ^ (A e - p /A e+ »)+e-^+^ 
P ~ (A e - p /A e+n )e"--^+« + l' 

where n e = u e /T is the electron degeneracy parameter 
and £ = (m„ - m p ) /T = Sm np /T w 1.293 MeV/T is the 
neutron-proton mass difference divided by temperature. 
Eq. 1)58(1 is approximate because it assumes identical neu- 
trino and plasma temperatures. We have also neglected 
the neutron decay/three-body capture processes of Eq. 
l|55jl. The expression in Eq. 1(58(1 is generally true for 
Fermi-Dirac leptonic energy distribution functions, even 
if the neutrinos and electrons/positrons are not in true 
thermal and chemical equilibrium. If and only if chemical 
equilibrium actually obtains (or did obtain at some early 
epoch) are we guaranteed to have fi e — /i„ e = jj, n — u p , 
where u n and u p are the neutron and proton total chem- 
ical potentials, respectively, and only in this case does 
Eq. JHHI) reduce to 

- « e (lMe-IJ.u e -Sm np )/T_ fags 
P 

With strict chemical equilibrium and with Fermi-Dirac 
energy distributions for all leptons, we could conclude 
from Eq. I(59|) . for example, that a positive chemical po- 
tential for electron neutrinos (i.e., an excess of v e over 
v e ) would suppress the steady state equilibrium neutron- 
to-proton ratio relative to that for n Uls = 0. This be- 
haviour follows also from a straightforward application 
of Le Chatlier's principle to the processes in Eqs. I(53|) 
and 1(54(1 . A decrease in the neutron abundance trans- 
lates, in turn, into a decrease in the predicted He yield. 

However, if the neutrino distribution functions are 
modified by active-sterile neutrino conversion, v a ^ v s , 
then the resulting active and sterile neutrino distribution 
functions would not be Fermi-Dirac in character and we 
could not employ Eq. I(59|) to determine the neutron-to- 
protion ratio in steady state equilibrium. Instead, we 
would be forced in this case to evaluate and follow the 
rates directly. 

We solve Eq. 1(56(1 numerically, assuming a homoge- 
neous and isotropic FLRW universe. In these calculations 
we take the co-moving entropy density to be conserved 



and thereby calculate self consistently the e densities 
and the neutrino and plasma (photon/e ± ) temperatures. 
All neutrino energy densities are handled correctly for all 
assumed lepton numbers and sterile neutrino populations 
and, therefore, the expansion rate is also self consistently 
calculated. At each time step in these computations we 
employ appropriate neutrino, antineutrino and e distri- 
bution functions and calculate the lepton capture rates 
Ajy c n, A e - p , Ap eP , A e +„, and the free neutron decay rate 
Andccay For convenience we adjust the Coulomb wave 
correction factor (see Appendix A) for the free neutron 
decay to be about (G) « 1.0227 to give a vacuum (un- 
blocked) neutron lifetime of 888 s. We employ a Coulomb 
wave correction (G) = 1 for all of our lepton capture rates 
and we adopt an effective ft- value for all weak reactions 
ft w 10 3 035 , so that there is roughly a ~ 2% incon- 
sistency in overall coupling between the lepton capture 
processes on the one hand and free neutron decay on the 
other in our calculations. This is another reason why our 
calculations of n/p can be used only to compare trends 
in various cases with and without spectral distortion and 
not as quantitative nucleosynthesis predictions. 

To gauge the effect of active-sterile neutrino flavor con- 
version on the evolution of the n/p ratio we employ a 
modified forced continuous resonance sweep scenario in 
either the v e — > v s or v e — > v s channels. We use Eq. I(25|) 
but fix the potential lepton number at half the initial 
value. This will leave a low energy "hole" with zero v e 
(or v e ) population, as in Fig. Q), that will grow with the 
expansion of the universe until e c . . is reached. The final 
high energy edge of the "hole," e c . . is determined from 
the initial potential lepton number as discussed in the 
previous sections. 

The modification of the lepton capture rates in the 
forced continuous resonance sweep scenario is discussed 
in Appendix A. Consider, e.g., v e — > v s . Because in this 
case there are now fewer u e 's, electron capture on protons 
will be less Fermi blocked and, hence, the capture rate, 
A e - p , will be larger. By the same token, fewer v e 's will 
translate into a reduction of the v e capture rate on neu- 
trons, \ v< , n - Note that both a larger value for A e - p and 
a smaller value for A t/c „ go in the direction of increasing 
the neutron-to-proton ratio in weak steady state equilib- 
rium. This is obvious from Eq. ((57|l . Likewise, this trend 
in the rates will similarly make itself felt in the general 
solution of Eq. ((56(1 and the net result will be an increase 
in n/p at T a . 

Of course, we argued above that the actual neutrino 
or antineutrino spectral distortions are likely to be quite 
different from the simple ones that continuous resonance 
sweep would give. Since we cannot calculate the actual 
neutrino spectra, we cannot give a quantitative calcula- 
tion of nucleosynthesis yields. We seek here to give rough 
guidelines as to where one might expect significant neu- 
trino spectral distortion modification of helium yields. 

The forced continuous resonance sweep picture will 
serve to get the general features of the rate effects. Fur- 
thermore, it sometimes will do this in a conservative 
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manner, especially for low lepton numbers where e c . . is 
small, i.e., e c o . < 3. As shown in Appendix A, the weak 
cross sections weight the square of the neutrino energy. 
A more realistic "picket fence" neutrino spectrum will 
remove population at higher neutrino energies than will 
the continuous sweep cut-off spectrum and will, there- 
fore, sometimes result in larger suppression of neutrino 
capture rates and increases in e ± capture rates. Like- 
wise, active-active mixing will partially fill in the "hole" 
in the neutrino spectra, but at the cost of pushing the 
neutrino spectral deficit (relative to a thermal spectrum) 
to higher energy. 

Let us consider a particular active-sterile neutrino con- 
version scenario in the channel v e — > v s . In this example 
we take each of the three neutrino flavors to have a lep- 
ton number near or at the maximum allowed without 
spectral distortion. We take L v = L„ x = 0.15, corre- 
sponding to degeneracy parameters n Uii = r\ Vr ~ 0.219, 
and take L Ve « 0.0343, corresponding to electron neu- 
trino degeneracy parameter r\ Vi! = 0.05. This will give an 
initial potential lepton number in the v & — > v s transfor- 
mation channel, 

^initial = + L ^ + L ^ ^ 3 g 8 (6Q) 

In this case the difference in neutrino energy density over 
the zero lepton case is only ~ 0.2%. Therefore, the ex- 
pansion rate of the universe at Weak Freeze Out in this 
case will differ from the standard BBN model by only 
~ 0.2%. Therefore, the expansion rate by itself would 
give a negligible difference in neutron-to-proton ratio be- 
tween the case with the lepton number in Eq. I|60|) and 
the zero lepton number, standard BBN case. 

Nevertheless, we calculate the n/p ratio completely self 
consistently as described above, beginning with steady 
state equilibrium at T = 3 MeV. The results are shown 
in Fig. 0. The lower solid curve is the n/p ratio in the 
case with the lepton numbers but with no active-sterile 
neutrino flavor transformation and, therefore, no sterile 
neutrino population and no active neutrino spectral dis- 
tortion. As expected, the n/p ratio is suppressed relative 
to a zero-lepton number standard BBN calculation owing 
to the large positive v e degeneracy parameter rj Vc = 0.05. 
At the lowest temperature on the plot, T ss 80keV (this 
is near T a in an actual nucleosynthesis calculation for 
baryon number n = 6 x 10~ 10 ) we have n/p w 0.133, 
corresponding to a rough 4 He yield X a ss 23.5%. A sim- 
ilar calculation but with zero lepton numbers (this is the 
standard cosmological model case) gives n/p w 0.141, 
corresponding to A° w 24.7%, so our lepton number- 
only case corresponds to about a 5% decrease in helium 
yield, as expected. 

The dashed line on Fig. is the equilibrium n/p ra- 
tio given by Eq. l(5^|) or Eq. l(5^|) (they give the same 
result in this case because all leptons have Fermi-Dirac 
distribution functions). We see that the lepton cap- 
ture rates become too slow to maintain equilibrium when 
T < 2 MeV, though the actual n/p ratio is still signifi- 
cantly influenced by these rates for temperatures greater 



than T ~ 300 keV. 

The upper solid curve in Fig. is the n/p ratio in the 
case where the same lepton numbers now drive v e — > v s 
in the modified forced continuous resonance sweep case 
as described above with Sni 2 = 10eV 2 . In this scenario 
the resonance will, eventually, sweep out to e c . . ~ 2.724. 
The upper dotted curve is the steady state equilibrium 
n/p ratio for this case, given by Eq. The "kink" 

in this curve (mirrored in a similar, smaller deviation in 
the solid curve) at about T w 1.2 MeV corresponds to 
the point where for this 5m 2 and (1/2) potential lepton 
number the resonance has swept far enough (i.e., near 
e c .o., see Eq. I15fl to significantly reduce the v e + n — > 
p + e~ rate. From then on the n/p ratio tracks higher 
than the lepton number-only case and at T w 80keV 
we have n/p 0.159, corresponding to X a rss 27.4%, an 
11% increase over the standard BBN zero-lepton number 
case, and a whopping nearly 17% increase over the lepton 
number-only case. (In fact, the bigger neutron number 
density in this case likely would lead to a slightly earlier 
assembly of alpha particles, i.e., a higher T a and, hence, 
an even slightly bigger 4 He yield.) 

In other words, the existence of a sterile neutrino that 
mixes with the v e completely altered the sign of the ef- 
fect of a net lepton number. The lepton number by itself 
would have given a comfortable reduction of the helium 
yield, whereas the v e spectral distortion in this case re- 
engineered this into an uncomfortable increase in the he- 
lium over the standard model. Despite the limitations 
of our calculations and approximations it is clear from 
this example that the existence of sterile neutrinos could 
alter the relationship between lepton number(s) and 4 He 
yield. 

Note that there would be significant alterations in 
the lepton number/helium-yield relation even if the res- 
onance did not sweep beyond e max . We have already 
argued that the resonance will find a way to sweep be- 
yond e max , given that the resonance condition can be met 
for non-continuous resonance sweep and that neutrino 
flavor evolution is likely quite adiabatic for the relevant 
conditions. However, for argument's sake, let us adopt 
the same example as above and with 8m 2 = 10 eV 2 , but 
now limit the resonance's progress to e max - For this ex- 
ample [L v = L Vt = 0.15 and L Ve s» 0.0343) we have 
£max ~ 1-461. The numerical Y e calculation in this case, 
using the same scenario as above, yields n/p ~ 0.138 at 
T « 80 keV, a nearly 4% increase over the lepton number- 
only case. 

Our spectral distortion effects can be dependent on 
Sm 2 , at least for the high lepton numbers adopted in the 
above example (L v — L Vt — 0.15 and L Ve « 0.0343). 
The middle solid curve in Fig. is the n/p ratio evolu- 
tion for this case with forced continuous resonance sweep 
out to e c . ., but now for Sm 2 — leV 2 . The lower dot- 
ted curve is the Eq. I|58|l steady state equilibrium n/p for 
this case. Since the resonance sweep is proportional to 
Sm 2 (see Eq. 1150 . we expect that the resonance sweep 
will not have progressed far enough to decrease the v e 
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capture rate significantly until a lower temperature than 
in the 8m 2 — 10 eV 2 case. That temperature is about 
T w 700 keV in this case. The result is that the rate 
reductions come in later here, where they are less effec- 
tive at lowering n/p, though there is still a hefty effect. 
At T 80keV we have n/p ~ 0.147, corresponding to 
X a w 25.6%, a nearly 4% increase over our standard 
BBN zero-lepton number case, and a 9% increase over 
the lepton number-only case. If we do the same calcula- 
tions but now for Sm 2 = 0.2 eV 2 we get X a « 24.7%, the 
same as the zero lepton number case and a 5% increase 
over the lepton number-only case, again a significant but 
smaller effect. 

However, we do not see this level of Sm 2 dependence in 
n/p alterations when the lepton numbers are small. This 
is because e ros ~ Sm 2 /C, so a low £ translates into more 
progress in resonance sweep for a given T and Sm 2 . 

A case in point is where L Ve = L v = L Vt = 0.01 
[y & degeneracy parameter rj Ve ~ 0.0146), corresponding 
to potential lepton number C e = 0.04 for the v e — > v s 
channel, with e c o . ss 0.96. This equilibrated case repre- 
sents a threshold: values of C e larger than this in fully 
equilibrated limit could produce significant (> 1%) in- 
creases in helium yield over the lepton number-only sce- 
nario, depending on the resonance sweep scenario and 
the efficacy of active-active transformation and residual 
neutrino down-scattering. The lepton number-only cal- 
culation for this case (as in Fig. 0) at T w 79 keV gives 
n/p w 0.139, or X a w 24.4%, a very slight decrease from 
the standard zero lepton number case discussed above. 
However, with 8m 2 > 0.2 eV 2 and v e ^ v s conversion en- 
abled, at T « 79keV we get n/p m 0.14, or X a m 24.6, a 
~ 1% increase over the lepton number-only case. Larger 
values of C e will give bigger discrepancies between the 
cases with and without sterile neutrino mixing in the 
forced continuous sweep limit, and may well do so in 
more realistic resonance sweep scenarios as well. 

Likewise, we can investigate the analogous limit for un- 
equilibrated cases (where active-active mixing is ineffec- 
tive) by employing n/p evolution calculations along the 
lines of those presented in Fig. J7J). We find that the case 
with L Ue = 0.001 and L v = L Ut = 0.01 (corresponding 
to C e = 0.022, e c . . « 0.76 and n Ve « 0.00146) is likely to 
give about a 1% increase in X a over the lepton number- 
only case in the forced continuous resonance sweep sce- 
nario. Again, depending on the actual resonance sweep 
history and the efficacy of neutrino scattering this repre- 
sents a warning point for the un-equilibrated cases. 

With these numerical calculations we are led to ask 
the following question: at what point are the positive 
lepton numbers big enough that when combined with a 
sterile neutrino and concomitant spectral distortions we 
exceed the classic observationally inferred helium limits 
discussed above (X^ m w 25%)? In the forced continu- 
ous resonance sweep scenario, and taking active-active 
mixing as efficient, we find that limit occurs for equili- 
brated lepton numbers near L„ e = L v = L„ T = 0.09 [v e 
degeneracy parameter r\ Ve ~ 0.131), corresponding to po- 
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FIG. 8: Constraints as in Fig. (|SJ and now with BBN "guide- 
lines," see text. The double-dot-dashed line gives an estimate 
of the limiting 4 He yield in a no active-active mixing scenario 
where L e = 0.001 and L M = L T = 0.01. The dashed line gives 
a threshold in potential lepton number for all cases where neu- 
trino lepton numbers are equilibrated, including those with ef- 
ficient active-active mixing. Beyond this threshold there could 
be significant alterations in the helium yield stemming from 
non-thermal neutrino distributions. The upper dash-dot line 
gives a rough estimate of the limiting potential lepton number 
for equilibrated cases with inefficient active-active mixing. 



tential lepton number C e = 0.36 for the v e —> v s channel, 
with e c o . w 2.57. The lepton number-only calculation 
for these parameters gives n/p 0.121 or X a w 21.6% 
at T w 79keV; whereas, active-sterile transformation 
in this case with 5m 2 — 10 eV 2 yields at this temper- 
ature n/p ~ 0.143 or X a w 25%. However, this drops to 
X a w 24.6% for Sm 2 = 5 eV 2 and is down to X a w 23% 
for 8m 2 — 0.2 eV 2 . For these lepton number parameters 
we note that 8m 2 > 5 eV 2 may already be ruled out by 
closure constraints. 

These rough cautionary guides as to where one might 
expect sterile neutrino-generated spectral distortions to 
become important are shown superposed on the closure 
fraction constraints in Fig. (JSJ) . Our calculations are only 
very general guides, as discussed above, as we cannot fol- 
low the resonance sweep physics in detail and, therefore, 
we cannot be quantitative about nucleosynthesis yields. 

One example of a significant deviation from the forced 
resonance sweep picture is where active-active neutrino 
mixing is efficient and neutrino lepton numbers equili- 
brate rapidly as resonance sweep occurs. In particular, 
in Cases 1 and 2 discussed in the last section, the neutrino 
populations in the "hole" are not zero, with the conse- 
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quence that for a given e c . Q . the rates are not as affected 
as in the cases discussed above. For exampie, for initial 
lepton numbers L Ue — L v = L Vt =0.1, corresponding 
to n Ve ps 0.146, the lepton number-only case steady state 
n/p ratio from Eq. at T = 0.7 MeV is n/p ps 0.136. 
By contrast, when we allow sterile neutrino production in 
Case 1, the resulting non-thermal v e spectrum increases 
the steady state value to n/p « 0.15. In Case 2 this goes 
to n/p ps 0.145 while in Case 3 it is only n/p ps 0.144. 
These examples serve to illustrate the change in lepton 
capture rates between Cases 1, 2, and 3. This gives a 
rough guide as to how spectral distortion effects on BBN 
might decrease in calculations with the full Y e evolution 
as we go from Case 1 toward Case 3. 

Likewise, the cases with the highest values of 5m 2 
coupled with lower values of £ in Fig. JHJ may already 
have experienced considerable resonant active-sterile con- 
version and concomitant lepton number depletion before 
Weak Decoupling, while active neutrino scattering down- 
scattering was still effective. This would either wash out 
much of the "hole" in the neutrino distribution functions 
or leave a hole of reduced energy width at higher energies. 

What about conversion of D e to steriles, i.e., D e — ► D s l 
This process can be matter-enhanced when the overall 
potential lepton number is negative, C e < 0. It will be 
exactly analogous to the positive potential lepton num- 
ber case, at least as far as the neutrino flavor conversion 
and the resonance sweep physics goes. This close anal- 
ogy ends, however, when it comes to the lepton capture 
reactions. 

If the resonance sweeps adiabatically in the forced reso- 
nance sweep scenario out to a scaled antineutrino energy 
e c .o., a hole will be left in the D e distribution, in com- 
plete analogy to the cases discussed above. Fewer i> e 's 
will translate into a decreased antineutrino capture rate, 
v e + P —> n + e + , and an increased rate of positron cap- 
ture, e + + n —>■ p + D e . These rate modifications both go 
in the direction of decreasing the neutron-to-proton ratio 
n/p, as is obvious from Eq. fBTjl . 

One might think at first that simply changing the sign 
of the potential lepton numbers given in the above ex- 
amples could result in a significantly suppressed n/p at 
T a . This is not necessarily correct however, because the 
threshold in the reaction v e + p — ► n + e + plays a cru- 
cial role. As can be seen in the rate integrals given in 
Appendix A, in this channel a v e must have an energy 
m excess of the threshold, E Pe > E^ TCsh to be captured. 
The threshold is El h J csh = + m ^2 _ pgQ4 MeV. 

In fact, in the forced continuous resonance sweep sce- 
nario, unless e c . . > E£ hresh /T, there will be no modifi- 
cations in the v e capture rates. Likewise for the inverse 
process of positron capture on neutrons, e + +n — > p + D e . 
In this reaction, there will be no alteration of the final 
state v e blocking factor unless e c . . > E^f esh /T. 

For example, consider the equilibrated case with L„ c = 
L Ufi = L Vt = -0.01 {rj Ve ps -0.0146). This gives the op- 
posite sign potential lepton number, C e ps —0.04, from 
the analogous positive lepton number case considered 
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FIG. 9: Constraints and guidelines as in Fig.s JSJ & (|5J, 
but now for negative values of potential lepton number. The 
dashed line gives a rough threshold on potential lepton num- 
ber in the equilibrated regime ( L e = L M = L T ) beyond 
which significant alterations (reductions) in 4 He yield stem- 
ming from Ve. spectral distortion might be expected. 



above. The forced continuous, adiabatic resonance sweep 
scenario would give for this case e c o . « 0.96. A calcula- 
tion of Y e with temperature as in Fig. JJJ) for the lepton 
number-only version of this case gives n/p ss 0.148 at 
T pa 79keV, or roughly X a ss 25.7%, an increase in he- 
lium yield over the standard zero lepton number case, as 
expected. However, when we now do the same calcula- 
tion but with v e — > v s conversion with Sm 2 > 0.2 eV 2 
we obtain at the same temperature n/p « 0.143, corre- 
sponding to X a ps 25%. This is a ~ 3% reduction over 
the lepton number-only case. 

This negative potential lepton number value again sig- 
nals a threshold: negative potential lepton numbers in 
the fully equilibrated limit larger in magnitude than this 
could give significant modification of the relationship be- 
tween lepton number and helium yield. This guideline is 
shown in Fig. © in the same manner as for the guides for 
positive potential lepton numbers. Again, the warning as 
to the very rough nature of our guidelines owing to uncer- 
tain resonance sweep physics and the efficacy of neutrino 
down-scattering and active-active mixing applies here as 
well as in Fig. ©. 

Because of the threshold issue in the negative potential 
lepton number cases it is not possible to find parameters 
that could provide a large suppression of 4 He over the 
zero lepton number case, at least in the fully equilibrated 
limit. Partly this is because large lepton number magni- 
tude would require a large dm 2 in order to have enough 
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resonance sweep early enough to affect the lepton cap- 
ture enough to drop n/p significantly. The large dm 2 at 
these large, negative C e values seem to be in conflict with 
closure constraints. 



V. DISCUSSION AND CONCLUSIONS 

The most general conclusions that can be drawn from 
this work is that the existence of one or more light sterile 
neutrinos could (1) alter the relationship between neu- 
trino chemical potential and primordial nucleosynthesis 
yields, and (2) leave both active and sterile neutrinos 
with nonthermal, distorted energy spectra. Attempts to 
constrain an "LSND" sterile neutrino based on conven- 
tional degenerate primordial nucleosynthesis considera- 
tions, as well as attempts to reconcile this neutrino with 
BBN limits via a primordial lepton number are now sus- 
pect. However, obtaining the detailed relationship be- 
tween lepton numbers, active-sterile neutrino mixing pa- 
rameters and light element nucleosynthesis yields that 
would be required to effect reliable constraints is beyond 
the scope of the work presented in this paper. 

This is mostly a consequence of another discovery 
made in this work: MSW resonances cannot sweep 
smoothly and continuously beyond e max . We showed 
that one way the resonance condition and adiabatic con- 
version criterion can be met beyond this point is if the 
resonance skips to higher energies a number of times un- 
til the initial lepton number is depleted. Partly for this 
reason resonance evolution is likely to be quite compli- 
cated. Complication will also arise because active-active 
neutrino mixing can be efficient and can occur simulta- 
neously with active-sterile transformation. 



In any case, neutrino flavor conversion is almost in- 
evitable once the resonance condition is met. This is 
on account of another insight presented in this paper: 
the highly adiabatic nature of neutrino flavor evolution 
through MSW resonances for the neutrino mass/mixing 
parameters of most interest and for the conditions in the 
post Weak Decoupling early universe. 

If the mini-BooNE experiment sees a positive signal, 
confirming the existence of light sterile neutrinos, we will 
be forced to confront the problems posed in this pa- 
per. Likewise, future progress in improving the precision 
and confidence in the observationally inferred primordial 
helium abundance coupled with CMB and large scale 
structure-derived limits on neutrino collisionless damp- 
ing scales could give us hints about active-sterile neutrino 
mixing physics in the early universe. 

Acknowledgments 

The work of G.M.F. was supported in part by NSF 
grant PHY-00-99499, the TSI collaboration's DOE Sci- 
DAC grant, and a UC/LANL CARE grant at UCSD; 
K.A. was supported by Los Alamos National Labora- 
tory (under DOE contract W-7405-ENG-36) and also ac- 
knowledges the UC/LANL CARE grant. N.B. was sup- 
ported by Fermilab, which is operated by URA under 
DOE contract No. DE-AC02-76CH03000, and was addi- 
tionally supported by NASA under NAG5-10842. This 
research was also supported in part by the NSF under 
Grant No. PHY99-07949. The authors thank the ECT* 
and the INT at the University of Washington for hospi- 
tality. We thank Ray Volkas and Robert Foot for discus- 
sions. 



[1] K. Eitel, New J. Phys. 2, 1 (2000). 

[2] G. McGregor [MiniBooNE Collaboration], AIP Conf. 

Proc. 655, 58 (2003). 
[3] H. Murayama and T. Yanagida, Phys. Lett. B 520, 

263 (2001); G. Barenboim, L. Borissov, J. Lykken and 

A. Y. Smirnov, JHEP 0210, 001 (2002). 
[4] M. C. Gonzalez-Garcia, M. Maltoni and T. Schwetz, 

Phys. Rev. D 68, 053007 (2003). 
[5] K. N. Abazajian, Astropart. Phys. 19, 303 (2003); P. Di 

Bari, Phys. Rev. D 65, 043509 (2002) [Addendum-ibid. 

D 67, 127301 (2003)]. 
[6] R. H. Cyburt, B. D. Fields, K. A. Olive and E. Skillman, 

Astropart. Phys. 23, 313 (2005). 
[7] K. A. Olive and E. D. Skillman, Astrophys. J. 617, 29 

(2004). 

[8] S. Hannestad, JCAP 0305, 004 (2003); . Pierce and 

H. Murayama, Phys. Lett. B 581, 218 (2004). 
[9] M. Tegmark et al. [SDSS Collaboration], Phys. Rev. D 
69, 103501 (2004). D. N. Spergcl et al, Astrophys. J. 
Suppl. 148, 175 (2003). 
[10] U. Seljak et al, Phys. Rev. D 71, 103515 (2005). 
[11] G. M. Fuller and D. N. Schramm, Phys. Rev. D 45, 2595 



(1992). 

[12] Z. Chacko, L. J. Hall, S. J. Oliver and M. Perelstein, 

Phys. Rev. Lett. 94, 111801 (2005). 
[13] Z. Chacko, L. J. Hall, T. Okui and S. J. Oliver, Phys. 

Rev. D 70, 085008 (2004); R. N. Mohapatra and S. Nasri, 

Phys. Rev. D 71, 053001 (2005). 
[14] J. F. Beacom, N. F. Bell and S. Dodelson, Phys. Rev. 

Lett. 93, 121302 (2004). 
[15] R. Foot and R. R. Volkas, Phys. Rev. Lett. 75, 4350 

(1995). 

[16] S.P. Mikheyev and AYu. Smirnov, Yad. Fiz. 42, 1441 
(1985) [Sov. J. Nucl. Phys. 42, 913 (1985)]; L. Wolfen- 
stein, Phys. Rev. D 17, 2369 (1978). 

[17] A. D. Dolgov, S. H. Hansen, S. Pastor, S. T. Petcov, 
G. G. Raffelt and D. V. Semikoz, Nucl. Phys. B 632, 363 
(2002); K. N. Abazajian, J. F. Beacom and N. F. Bell, 
Phys. Rev. D 66, 013008 (2002); Y. Y. Y. Wong, Phys. 
Rev. D 66, 025015 (2002). 

[18] V. Barger, J. P. Kneller, P. Langacker, D. Marfatia and 
G. Steigman, Phys. Lett. B 569, 123 (2003); 

[19] J. P. Kneller, R. J. Scherrer, G. Steigman and 
T. P. Walker, Phys. Rev. D 64, 123506 (2001). 



21 



[20] D. P. Kirilova, Astropart. Phys. 19, 409 (2003). 

[2lj K. Abazajian, X. Shi and G. M. Fuller, 

arXiv:astro-ph/9909320 
[22] M.Yu.Khlopov and S.T.Petcov Phys. Lett., B99, 117 

(1981); Ya.B.Zeldovich and M.Yu.Khlopov Usp. Fiz. 

Nauk 135, 45 (1981). [English translation: Sov. Phys. 

Uspekhi 24, 755 (1981).] 
[23] M. Cirelli, G. Marandella, A. Strumia and F. Vissan i, 

Nucl. Phys. B 708, 215 (2005) |arXiv:hep -ph/0403158 . 
[24] D. Notzold and G. Raffelt, Nucl. Phys. B 307, 924 (1988). 
[25] X. Shi and G. M. Fuller, Phys. Rev. Lett. 82, 2832 (1999). 
[26] Y.-Z. Qian and G. M. Fuller, Phys. Rev. D 51, 1479 

(1995) . 

[27] A. B. Balantekin and G. M. Fuller, Phys. Lett. B471, 
195 (1999). 

[28] G. C. McLaughlin, J. M. Fetter, A. B. Balantekin and 

G. M. Fuller, Phys. Rev. C 59, 2873 (1999). 
[29] M. Sorel, J. M. Conrad and M. Shaevitz, Phys. Rev. D 

70, 073004 (2004). 
[30] J. M. O'Meara, D. Tytler, D. Kirkman, N. Suzuki, 

J. X. Prochaska, D. Lubin and A. M. Wolfe, Astrophys. 

J. 552, 718 (2001). 
[31] K. A. Olive, E. Skillman and G. Steigman, Astrophys. J. 

483, 788 (1997). 
[32] Y. I. Izotov and T. X. Thuan, Astrophys. J. 602, 200 

(2004). 

[33] Y. Z. Qian, G. M. Fuller, G. J. Mathews, R. Mayle, 
J. R. Wilson and S. E. Woosley, Phys. Rev. Lett. 71, 
1965 (1993). 

[34] M. S. Smith, L. H. Kawano and R. A. Malaney, Astro- 
phys. J. Suppl. 85, 219 (1993). 
[35] C. Y. Cardall and G. M. Fuller, Astrophys. J. 472, 435 

(1996) . 

[36] G. M. Fuller, W. A. Fowler, M. J. Newman, Astrophys. 
J. 293, 1 (1985). 



APPENDIX A: WEAK RATES WITH 
NON-THERMAL NEUTRINO ENERGY 
SPECTRA 

In this Appendix we calculate the forward and reverse 
rates of the processes in Eqs. 11531 & (|541 for the cases 



Here the integration variable in both equations is the 
scaled v e energy, x = E v<t jT. The final state electron 



where, respectively, all neutrinos below energy E Ve = Te 
or antineutrinos below energy Ep e = Te are converted 
to sterile species. We provide estimates of these rates 
in terms of standard relativistic Fermi integrals. We 
also discuss how these rates would be modified if the 
MSW resonance does not sweep smoothly and continu- 
ously (and adiabatically) through the low energy neu- 
trino or antineutrino distribution function, but instead 
skips to higher energies. The rate modifications for Cases 
1 and 2 in the efficient active-active neutrino mixing limit 
will be different, of course, because in those scenarios 
the "holes" in the neutrino distribution functions are not 
empty. Though the rate formulae presented here are not 
valid for these cases, they still give a general idea of how 
the lepton capture/decay rates depend on spectral dis- 
tortion and thresholds. 



If there is no active-sterile conversion and all neutrino, 
nucleon, and charged lepton distribution functions are 
thermal in character, the v e capture rate on neutrons 
is A° e „. By contrast, we will denote as X Ven the actual 
electron neutrino capture rate when the same thermody- 
namic conditions obtain, but now where v e J s have been 
converted to sterile species up to scaled energy e as out- 
lined above. If all neutrino, nucleon, and charged lep- 
ton energy distribution functions are at least piece-wise 
Fermi-Dirac or zero, these rates can be written [3(|, re- 
spectively, as 



^jdx, (Al) 
- J dx (A2) 



energy is E e = T (x + £ np ). The other notation in these 
expressions is as defined above and £ np = Q np /T with 
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Qnp — 5m np . There is no threshold for v e energy in this 
reaction channel. The temperature and matrix clement- 
dependent factor in both rate expressions is 



A 



{ft) \m e c z 



f 1.835 x lO^s- 1 



(A3) 



where {ft) is the effective ft- value as defined in Ref. |3|j 
and is roughly log 10 ft w 3.035 for free nucleons, while 
(G) is the average Coulomb wave correction factor (also 
defined in Ref. |2g) with G = F(Z, E e )E e /p e and where 
F(Z,E e ) is the usual Fermi function for nuclear charge 
Z and final state electron energy E e . For the relativistic 
leptons considered here (the lowest electron energy is rs 
Q np m 1.3 MeV), (G) « 1, though we note that (G) in 
the no-transformation case is slightly larger than that for 
the case with the e cut-off on account of the lower energy 
electrons present in the phase space integral in the former 
case. (Electrons are "pulled in" to the proton because of 
Coulomb attraction, making for a larger overlap.) 

The second approximation in Eq. l|A2fl gives A„ e „ as a 
sum of differences of relativistic Fermi integrals. In this 
expression the effective v e and e~ degeneracy parameters 
are defined as ?y° ff = T) Ve — e and 7y° ff = r] e — S, respectively, 
where 6 = e + £ np . Also in Eq. (|A2|) we define 04 
1 . while «3 = 2{e + S), and «2 = (e + <5) 2 + 2eS, with 
oil = 2e(5 (e + S) and «o = e 2 <5 2 . Note that as e — > 0, 
both expressions in Eq. 1A2I) approach A°_ in Eq. 1A1JI . 

e P 

It is obvious that for nonzero e the v e capture rate on 
neutrons will be reduced over its no-transformation value, 
Ay c n < A° n . The rate A„ e „ is illustrated in Fig. 02] as a 
function of e or L Ue . 

The rate for the corresponding reverse process of elec- 
tron capture on protons, e~ + p — > n + v e , will be in- 
creased if some VeS are transformed to sterile states, as 
there will be less final state v e blocking in this case. For 
a Fermi-Dirac distribution of electrons, and in terms of 
an integral over electron energy E e , this rate is 



(G) In 2 

(/*)(m e c 2 ; 



El {E e - Q 



np J 



,E e /T- V e + 1 



[l-S v ,]dE e , 
(A4) 

where S Ve is the energy-dependent v e occupation proba- 
bility, 



s Ve = 



1 



E„ e /T- Vue + 1 



for E V JT < e, (A5) 
for E U JT> e. (A6) 



Here the v e energy is E„ e — E e — Q np on account of the 
threshold, Q np . 

It is convenient to re-write the rate in Eq. I|A4() as 
an integration over neutrino energy scaled by tempera- 
ture, x = E U JT, and as a sum of contributions from low 
neutrino energy with no final state blocking, and higher 




FIG. 10: Rate \ Vt , n in s _1 for the process v E + n — > p + e~ at 
temperature T — 0.7 MeV as a function of e and/or L Ve in the 
smooth and continuous resonance sweep limit and for the case 
of complete active neutrino equalization {L Ve — L Vfl — L Vt ). 
The solid curve gives the rate for no sterile neutrino con- 
version, thermal v e distribution, but with the u e chemical 
potential appropriate for the corresponding e value. The dot- 
dashed curve gives the rate with active-sterile neutrino con- 
version and corresponding non-thermal character for the u e 
energy distribution function. 



final state neutrino energy where there is non-zero Fermi 
blocking, 



A, 



\ low 

A e~p 



A 



high 



(A7) 



The first of these rate contributions can be approximated 
by 



low 



A 



(■£ F £np) 



1 



-dx. 



(A8) 



Just as for v e capture, the average Coulomb wave correc- 
tion factor will be lower (closer to unity) with increasing 
e. Again this has to do with the enhancement of the low 
energy electron probability density near the proton. As 
above, we can represent the rate contribution in Eq. 
in terms of standard relativistic Fermi integrals, 



A 1 ,™ « A [F 4 ( Ve - U) + ^n P F 3 (r? e - e„ P ) + £ P F 3 ( Ve - 



A^/3, 

n=0 



(A9) 
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FIG. 11: Rate \ ep in s _1 for the process e~ +p —* n + v e at 
temperature T — 0.7 MeV as a function of e and/or L u<s in the 
smooth and continuous resonance sweep limit and for the case 
of complete active neutrino equalization (L Ve = L Vfl — L Vt ). 
The solid curve gives the rate for no sterile neutrino con- 
version, thermal v e distribution, but with the v e chemical 
potential appropriate for the corresponding e value. The dot- 
dashed curve gives the rate with active-sterile neutrino con- 
version and corresponding non-thermal character of the v R 
energy distribution function. 



Ae- S p « A[l 

« A[l 

where the notation is as above and where the fa are as 
defined above for Eq. (|A9|I . In summary, a hole in the low 
energy v e distribution results in a lower value for X Uan , a 
higher value for A e - p and, hence, an increased n/p ratio. 
The rate \ e - p is illustrated in Fig. ^] as a function of e 
or L Ue . 

By contrast, conversion up to scaled energy e = Ep e /T 
of iVs to sterile neutrinos, T> e — > D s , would result in a 
lower value of the neutron-to-proton ratio and, hence, 



where fa = 1, and where fa = 2(e + <5), while fa 
i + Sf + 2eS and fa = 2e8(e + 6), with fa = e 2 8 2 . Here 
we define 5 = e + £ np . 



The physical interpretation of this expression for \ °™ p 
is clear if it is recalled that the v e energy is E Ve = 
E e — Qnp, implying that the "effective final state neutrino 
degeneracy parameter" is rj e — £ np for the no-conversion 
case, and ry e — £„ p — e with conversion of z/ e 's to steriles. 
Of course, as e — > 0, the rate contribution from the (final 
state v e ) unblocked portion of the phase space approaches 
zero, A low — > 0. The second of the rate contributions in 

' e p 

Eq. 1A7|I can be approximated as 



(A10) 



a lower 4 He yield. This is because a low energy deficit 
in the D e distribution would lead to a decreased rate for 
v e + p — > n + e + and, on account of less blocking, an in- 
creased rate for the reverse process. Handling the energy 
threshold for these reactions is, however, somewhat more 
complicated than for v e and e~ capture. 

Using much the same notation as above, we can ap- 
proximate the rate for v e + p — > n + e + as 



I 

" 1 jf + 6.)* ( e . +( J-„. + 1 - ^Tl) * 

4 

" ^ Pn [F n (Ve ~ Cnp ~ e) - F n - e)], 

n=0 

I 



J 



« A[l - e*"-^-**]- 1 J2 (vf) ~ F n (rjf)] . 

I 

The integration variable in the first of these equations is x — E^^/T, and the final state positron energy will be 



24 



E e + = T (x — £ np ). The scaled energy threshold in these 
expressions is 

7thrcsh = €n P + m e for £ np + m e > e (A12) 
7throsh = e for e > £ np + m e 

where m e = m e c 2 /T. It is clear that transformation of 
i> e 's with energies below the threshold energy Q np + m e c 2 
does not affect the rate. In the second approximation in 
Eq. (|AII() , the effective D e degeneracy parameter is rjf 1 = 



T\y B — e, while the effective positron degeneracy parameter 
is ry| ff = t; np — r/ e — e. (Since electromagnetic equilibrium 
always obtains here, the positron and electron degeneracy 
parameters have equal magnitudes and opposite signs, 
rj e + = —rj e .) If we define a = 2e — £ np and b = e (e — £ np ), 
then the coefficients a n are: 0:4 — 1; S3 = 2a; <52 = 
a 2 + 26; ai = 2a6; and ao — b 2 . 

Utilizing the same quantities and notation as in Eq. 
(jAllfl . the rate for the reverse process of positron cap- 
ture, e + + n — ► p + D e , can be written as 



A 

1 _ e Ve-(np+Vn 



X {x £np) 



1 



1 



\dx + A 



aX-£, np + 7] e _|_ I 



dx. (A13) 



Again we see that if e < £„ p + m e , then from Eq. (|AI2(I 
the threshold is 7thrcsh = Cnp + m e and the neutrino flavor 
conversion will have no affect on the rate. In this case, 
the second term of Eq. i|AI3(l will vanish and the first 
term will be the rate with no neutrino conversion. The 
full rate expression in Eq. (|A13|) can be broken up into 
three parts, 

K+n — Af+ S * + A™ d „ + Ag+ rd , (A14) 

each of which can be rendered in terms of standard rela- 
tivistic Fermi integrals. 

Here A fi + S ^ corresponds to the first integral in Eq. 
(|AI3|) . It can be reduced to 

« i _ eV X P+VBs E a n [F n ( V f) - F n ( v f )] , 

n=0 

(AI5) 

where the a n are as defined for Eq. i|AII|) . the effective 
positron degeneracy parameter is ry| ff = — r\ e + £ np — e, 
and the effective v e degeneracy parameter in this case is 

off - 

til = Ve c - e- 



Note that the second integral in Eq. (|AI3|) is the sum 
Ae+^i + A*+ rd • The last term in this sum can be approxi- 
mated as 

4 

A*+ rd w a n F n {Up ~Ve- e) , (A16) 

n=0 

where the a n are the same as defined above for Eqs. (|A11|) 
& Eq. (|A1 5fl . In similar fashion we can express A s + d in 
terms of standard relativistic Fermi integrals, 

4 

A£ d „ « A E PnFn (-m e - Ve)- (A17) 

We define x = 2m e 4- £„ p and y — m e (m e + £ np ), with 
m e = m e c 2 /T. With these definitions we can write the 
j3 n in Eq. l{A"T7|) _as: (3 4 = I, while (3 3 = 2x, (3 2 = x 2 + 2y, 
fii = 2xy, and /3 = y 2 . 



